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Abstract 

Topics in the description of the properties of charmonium states are reviewed with an emphasis 
on specific theoretical ideas and methods of relating those properties to the underlying theory of 
Quantum Chromodynamics. 
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1 Introduction 



Time is carrying us farther away from that day of November 11, 1974, when the news of an unusual 
resonance found simultaneously at BNL[1J and at SLAC[2] has swiftly spread through high-energy 
laboratories around the globe. Well before the age of widespread instantaneous satellite communications 
and the Internet, the new resonance became known in Moscow within the same hour as it was publicly 
announced in Palo Alto. After the initial excitement, confusion and revelations, it became clear that 
the new J/ip resonance was the first to have been observed state of a system containing previously 
unknown (but anticipated) charmed quark and its antiquark: cc0. The new system, charmonium, in a 
close analogy with positronium or even with a hydrogen atom, was expected to contain a spectrum of 
resonances, corresponding to various excitations of the heavy quark pair. However, unlike its analogs 
governed mainly by the electrostatic Coulomb force, the properties of charmonium are determined by 
the strong interaction, so that the newly found system was, in a way, the simplest object for a study 
of the strong interactions. It was strongly hoped [1] that charmonium could play the same role for 
understanding hadronic dynamics as the hydrogen atom played in understanding the atomic physics. 
In a way, this has indeed been the case and the development of many methods in QCD is directly 
related to analyses of the properties of charmonium and of its heavier sibling bottomonium. 

Recently the physics of charmonium regained a great renewed interest due to the massive dedicated 
investigation by BES and CLEO-c and the studies using decays of B mesons and the radiative return 
technique at the B factories with a higher initial energy of the electron and positron beams. After a 
'dry speir of more than two decades during which no new states of charmonium have been found with 
any certainty, new observations discover charmonium and charmonium-related resonances at a rate that 
outpaces the ability of the theory to fit their properties in a consistent scheme. Furthermore, the data 
very strongly suggest that among the new resonances there are exotic four-quark states, possibly hybrid 
states with gluonic degrees of freedom in addition to the cc pair, and also loosely bound states of heavy 
hadrons - charmonium molecules. Thus it looks like that charmonium not only has provided us with a 
'hadronic atomic physics' but quite possibly also with a 'hadronic chemistry', and in its mature age of 
33 charmonium still offers us new intriguing puzzles. 

In this review some properties of the old and new states of charmonium and QCD-based methods 
of study of these properties are discussed. An emphasis is made on selected theoretical methods rather 
than on presenting the whole field and reviewing the data. An all-inclusive presentation can be found in 
the much larger review and an excellent update on the most recent data and the related theoretical 
developments is given in Ref.[B]. 

The topics in the spectroscopy of the traditional charmonium states are provided in Section 2, 
and the annihilation of and the radiative transitions between these states are discussed in respectively 
Sections 3 and 4. The Section 5 is devoted to hadronic transitions between charmonium levels and the 
related topic of the interaction of slow charmonium with hadronic matter. Finally, in Section 6 are 
discussed peculiar properties of the resonances with masses above the open charm threshold. 

2 The Spectrum of Charmonium States 

2.1 General Considerations 

The diagram of the known charmonium and apparently charmonium-related states is shown in Fig{Tl 
The quantum numbers and basic properties of most of the states in the charmonium family can be 
described within a simple picture of a nonrelativistic quark - antiquark pair cc. In this picture the 
states are characterized by the orbital angular momentum L, the total spin S of the quark pair, and 

^The second charmonium resonance was found just ten days after the [3]. 
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Figure 1: The known charmonium and charmonium- related resonances and some transitions 
between them. Also are shown (dotted lines) the thresholds for various pairs of charmed 
mesons. 



the total angular momentum J, which defines the spin of the state viewed as a particle. As usual, the 
total angular momentum is given by the vector sum of the orbital and the spin momenta: J — L + S. 

— * 

Likewise, the total spin S is determined by the vector sum of the quark and antiquark spins: S — Sc+Sc- 
Clearly, S takes the values and 1, thus splitting the four possible spin states of the pair into a singlet 
and a triplet. Furthermore, the excitation of the radial motion of the cc pair results in a spectrum 
of levels with the same L, S and J, and differing by the "radial" quantum excitation number rir with 
rir — corresponding to the lowest state in this spectrum. It is therefore customary to encode the values 
of these quantum numbers for each state of charmonium in the form of the symbol (rir + l^'^^'^^^Lj. 
The combination 25* + 1 conveniently indicates the spin multiplicity, while following the tradition from 
atomic physics the values of L, L = 0, 1, 2, 3, . . . are written as S, P, D.F, . . .. In this notation the lowest 
state with L = 0, S = and (necessarily) J = is represented as 1^5*0 {rjc resonance) while the first 
excited state with the same quantum numbers is 2^So {rj'J. 

The value of L determines the parity (P) for each of the states: P = (—1)''"+^, while L and S com- 
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bined also determine the charge conjugation parity: C = (—1)^+"^. Therefore the previously mentioned 
^5*0 states have quantum numbers J^'^' = O""*", while e.g. ^Si states have the same quantum numbers 
J^*" = 1 as the electromagnetic current, so that these states {J/ip, i/j', ...) can be produced as 
resonances in e^e~ annihilation. 

This briefly described standard nomenclature of the quark-antiquark states originates in the strictly 
nonrelativistic mechanics. Relativistic effects in the dynamics of quark and antiquark interacting with 
each other preserve most, but not all of it. Indeed, the conservation of the total angular momentum 
ensures that the states have definite J. On the contrary, the value of L is generally not preserved by 
the interaction. In particular, the operator of the so-called tensor forces S {S + 1) — 3{S ■ r){S ■ r)/r'^ 
does not commute with for the states with 5* = 1, so that a mixing of states with different L takes 
place. However the parity conservation requires that only the states with the same parity of L can get 
mixed, and the conservation of J then implies that the values of L for the mixed states can differ by at 
most two units. For instance a ^5*1 state can receive, due to relativistic effects and admixture of a ^Di 
state ('5 — D mixing'). It can be further noticed that the mixing in L is necessarily absent for certain 
states. Indeed all states with J = L are pure in L, both the spin-singlets, S = 0, and the spin-triplets, 
S" = 1, as are the ^Pq resonances. 

The applicability, in a certain extent, of a nonrelativistic description to charmonium has always been 
a source of interest to this system with the hope that in it the quark dynamics can be studied being 
not overly complicated by the relativistic effects. The significance of such effects in charmonium can 
be very approximately estimated already from the masses of the resonances, e.g. the mass difference 
AM between the ground I'^Si state {J/ip resonance) and its first radial excitation 2^5*1 {tp') in units of 
either of the masses provides an estimate of the relativistic parameter f^/c^: 



Such moderate, but not very small magnitude of the relativistic effects in charmonium in fact makes 
some of these effects visible in experiments and thus further expands the range of dynamical details 
that can be studied in the charmonium system. 

2.2 Potential Models 

2.2.1 Leading Nonrelativistic Approximation 

One widely used approach to describing charmonium is to consider its dynamics in analogy with atomic 
systems or positronium and to treat it in the nonrelativistic limit by means of a Schrodinger equation 
with a potential V{r) depending on the distance r between the quark and the antiquark. The relativistic 
effects up to the order v^/c^ can then be considered as perturbation due to relativistic terms in the 
potential as well as in the kinetic energy. The shape of the potential at short distances is determined by 
the perturbation theory in QCD. In the lowest order the exchange of (Coulomb) gluons between slow 
quarks is fully analogous to interaction in QED, so that for a color-singlet quark pair the interaction 
takes the Coulomb-like form: 



where is the QCD coupling. Once the scale dependence of this coupling is taken into account, the 
constant in Eq. ([2]) has to be replaced by the running coupling as(r). At distances longer than the 
charmed quark Compton wave length l/rric, the one- loop running is described by 



AM 




M 



0.2 . 






(3) 



5 



In higher orders of the QCD loop expansion the precise relation of this 'Coulomb-like' coupling to the 
running constant defined in a specific renormalization scheme, such as e.g. the MS scheme, is a matter 
of calculations, which have been carried to the two-loop level [71 |B] with some partial results [H] in three 
loops. 

The details of these fine calculations however are not of an immediate significance for charmonium. 
The reason is that the perturbative QCD is applicable only at short distances, which are much shorter 
than the typical spatial size of the charmonium states. At the relevant intermediate and long distances 
one has to resort to models for the interaction between quarks. Some guidance in constructing such 
models is provided by the general idea of quark confinement, which can be mimicked by a potential rising 
at long distances. The most popular choice of the confining behavior is a linearly growing potential: 
V{r) = ar. Such behavior originates in the idea of the contraction of the chromoelectric field between 
the quarks into a flux tube, giving a string-like binding. 

The interaction potential can also be studied as the energy of a static infinitely heavy quark - 
antiquark pair separated by the distance r. This quantity can be evaluated in terms of the Wilson 
loop[TO] by lattice QCD calculations. The numerical results of such analysis[Tl] produce a dependence 
of the static energy on r, which is in agreement with the Coulomb-like behavior at short distances and 
an approximately linearly rising potential at larger r. 

It should be noted however that unlike in QED, a potential approach to heavy quarkonium in QCD 
is formally justified only in the limit of very high quark mass: in tens to hundreds GeV. In this limit 
the low-lying bound states in the short- distance potential ([2]) are localized at short distances, where 
the perturbative potential description is applicable and the whole approach is thus self consistent. Once 
nonperturbative effects in QCD are taken into account such consistency becomes questionable. Already 
the leading corrections [121 HSl [E] to energies of the quarkonium levels in the limit of very heavy quarks 
do not correspond to any potential between the quark and the antiquark. The reason for such behavior 
can be readily understood [12]. Indeed, a potential implies an instantaneous interaction. Any non- 
locality of the interaction in time would contain characteristic time scales, that can be interpreted as 
the evolution time scales for additional degrees of freedom. Once such additional degrees of freedom 
come into play the system (quarkonium) can no longer be described by a potential, neither it can be 
considered at all as a two-body system. In reality an interaction between the quark and antiquark 
through a gluon field should necessarily invoke nonperturbative light degrees of freedom in QCD whose 
typical evolution scale is determined by the QCD infrared parameter Aqcd- The interaction through 
exchange of such field can thus be viewed as instantaneous inasmuch as the quark and antiquark are 
slow in this scale, i.e. the characteristic time of evolution of the quarkonium wave function is long in this 
scale. For charmonium one can estimate the time of evolution as an inverse of the typical energy spacing 
between the levels, e.g. M{il)') — M{J/iIj) ^ 590 MeV, which by any measure is certainly comparable 
with the QCD scale. 

In other words, there in fact is no parameter that would justify a QCD-derived description of 
charmonium or bottomonium as a two-body quark-antiquark system interacting through a potential. 
However due to some numerical reasons, which are yet to be understood, such simple picture works 
reasonably well, especially if it is not pushed to requirements of high accuracy, or to highly excited 
states of quarkonium. It is almost so that any smooth potential whose behavior resembles Coulomb 
at short distances and an approximately linear rise at large r reasonably well describes the properties 
of the observed charmonium resonances, after the parameters of the model are appropriately adjusted. 
Some of the models for the potential discussed in the hterature can be found in Refs. [TJl [ISl HTj . 

One of the most developed is the Cornell model[ISl dSl [H] which builds upon the simplest potential 
being just a sum of the Coulomb and linear parts. 



and adding then finer effects, such as the relativistic terms, resulting in the hyperfine and fine structures 
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of charmonium levels, and also, importantly, including the coupling of the cc system to pairs of charmed 
mesons, such as DD. The latter coupling effectively accounts for the fact that above the open charm 
threshold the charmed quarks and antiquarks do emerge inside the charmed hadrons, which effect is 
totally ignored in a pure potential model with confining interaction. 

Clearly, a leading nonrelativistic treatment can only describe gross features of the charmonium levels, 
i.e. without resolving the fine splitting between the states with the same L and S and different J and 
the hyperfine splitting between the spin-triplet and spin-singlet states. It can be noted that even at 
such approximate level of detail the resulting sequencing of the energies of the levels with different rir 
and L can provide useful constraints on the properties of the potential \^(r) j20]. 



2.2.2 Spin- dependent Forces 



The potential description extended to spin-dependent interactions results in three types of interaction 
terms that are to be added to the discussed leading nonrelativistic interaction: 



Vi{r) = VLs{r) {L-S)+VT{r) 



S{S+l) 



+ Vss{r) 



SiS + 1)-- 



(5) 



The spin-orbit, Vls: and the tensor, Vr, 
spin term, Vss, proportional to 2(s„ ■ 



terms describe the fine structure of the states, while the spin- 
) = 5'(5' + 1) — 3/2 gives the spin-singlet - triplet splittings. 



The interaction in Eq.([5]) arises among the v^/c^ effects in the nonrelativistic expansion and it generally 
requires additional mo del- dependent assumptions about the structure of the interquark forces. Within 
the phenomenological approach it is usually assumed that parts of the static potential (similar to 
that in Eq.Q) correspond to definite Lorentz structures of the relativistic interaction between the 
quarks [2T1 [221 H]. In other words, those structures correspond to an "exchange of something" with a 
definite spin between the quark and the antiquark. Then the short-distance Coulomb-like part of the 
static potential is naturally generalized as a limit of a vector type exchange: 

with u and v being the Dirac spinors for the quark and the antiquark, while the confining part has 
been treated in the literature as a part of a vector exchange |2H 1^ . or a scalar [51 [B], or a mixture of 
these|23j. With this choice of options restricted to a combination of only vector and scalar exchange, 
the spin-dependent terms in Eq.Q can be written in terms of the vector, Vv{r), and scalar, Vs{r), 
parts of the static potential by the standard Breit-Fermi expansion to order w^/c^[24j: 



Vj 
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2m^r 



dr 



dr 



(6) 



where m is the charmed quark mass. 
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r dr 



(7) 



AVv 



with A = being the three-dimensional Laplacian. 

It should be emphasized that by its nature the interaction in Eq.([5]) is a part of the v'^/c^ term in the 
nonrelativistic expansion. As such it can be used only in the first order, and by no means one should 
iterate this potential, or use it for an input in the Schrodinger equation. Therefore any accuracy of the 
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results found with potential of this type is intrinsically limited. Furthermore, the formulas ([6]) - ([8]) 
are only approximate even in perturbative QCD at short distances. In particular, the tree-level QCD 
potential ([2]), results, according to Eq.Q, in a point-like spin-spin interaction: 

which is the correct tree-level expression for the corresponding interaction in QCD. If one then improves 
the static potential by using the running coupling from Eq.(j3]), the formula ([8]) would produce terms 
extending to finite r and behaving as a^/r^. A real calculation [25] of the one- loop correction to the 
spin-spin interaction, however produces no such terms, and predicts that the hyperfine splitting in 
perturbative QCD is still proportional to the square of the wave function at the origin, |-?/'(0)p. 

Clearly, the point-like behavior of the spin-spin interaction can generally be invalidated in higher 
orders in perturbative QCD and also by nonperturbative dynamics. Moreover, the leading nonper- 
turbative effects in the limit of asymptotically heavy quarkonium|26] are not reduced to a point-like 
spin-spin interaction. Nevertheless, inspite of these reservations, the actual hyperfine splitting in char- 
monium closely resembles that produced by a short-distance interaction. Namely, the proportionality 
of the hyperfine splitting to |^/'(0)p implies that the mass gap between the ^5*1 and ^Sq states should be 
proportional to the e+e~ decay width of the vector ^5*1 resonance, while the hyperfine splitting in the P 
wave should be extremely small because of vanishing wave function at the origin. In phenomenological 
terms this implies the following relations: 

M{J/i;) - M(r/e) ^ TM/^) ' 

and 

M{K) ^ Mixcj) . (11) 

where M{xcj) = [5M(xc2) + 3M(xci) + M{xc2)]/Q is the 'center of gravity' of the ^Pj states which is 
not shifted by either the spin-orbital or the tensor interactions from Eq.([5]). According to the Tables [27] 
the ratio of the e~^e~ decay rates in the r.h.s. of Eq. ffTOl) is 0.45 ± 0.02, while the ratio of the mass 
splittings in the l.h.s. is 0.44 ± 0.04. Furthermore, the center of gravity of the spin-triplet Xcj states is 
at M(xcj) = 3525.36 ± 0.06 MeV. The spin-singlet ^Pi state, the he, has been sighted [28] as a resonance 
in the pp annihilation with the mass 3526.28 ± 0.18 ± 0.19 MeV, then in the same processpQ] at the 
mass 3525.8 ± 0.2 ± 0.2 MeV, and eventually [30] in the decays ^' ^ 7r° K at 3524 ± 0.6 ± 0.4 MeV, with 
the most recent improvement in the precision yielding[3T] the he mass of 3525.35 ± 0.19 ± 0.15 MeV. 
The data with the smallest claimed experimental errors point at an extremely small violation, if any, of 
the relation ([II]): M{hc) - M{Xcj) = -0.05 ± 0.19 ± 0.16 MeVfST]. Thus the simple relations 1^ and 
( ITTl) both hold amazingly well. In fact it is still a challenge for experiment to measure the violation of 
these relations, and even a greater challenge to correctly predict such violation theoretically. 

The spin-orbit and tensor terms in Eq.Q produce the fine structure of the charmonium levels, and 
the tensor term also gives rise to mixing of states with L differing by two units, such as ^5*1 — ^Di 
mixing. The phenomenological effects of the mixing are somewhat more subtle and will be discussed 
further in Sec. 3. 1.2. Here we concentrate on the fine structure of the ^Pj states Xcj- The shifts of the 
masses of the states with different J with respect to M(xcj) are given in terms of the averages < Vls > 
and < Vt > over the P wave coordinate wave function as 

6MCPo) = -2{Vls) + 2{Vt) , 6MCP,) = -{Vls) - (Vt) , SM^P^) = (Vls) + ^ (Vt) . (12) 

5 



Using the measured differences between the masses of the XcJ charmonium resonances [27] one can find 
the average values over the IP charmonium: 

{Vls) = ^ [5 M{xc2) - 3 M(xci) - 2 M(xco)] ~ 35 MeV , 
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{Vt) = ^ [M(xc2) - 3 M(xei) + 2 M(xco)] ~ -20 MeV . (13) 

This estimate illustrates that in order to describe the observed mass splitting between the Xcj states 
both the LS and tensor interactions are required with comparable strength. Furthermore, it can be 
noted that these forces, generated by a pure tree-level gluon exchange, as can be found from using 
the Coulomb-like potential Vy in the formulas (EI) and ([7]), have correct signs, in agreement with the 
estimate (fT3i) . but with a relative strength of the LS interaction enhanced in comparison with this 
estimate: Vis /Vt\couI = —3. Thus a certain reduction in the spin-orbit term due to a contribution of 
the Lorentz scalar potential Vs in Eq.([6]) is indeed helpful from this point of view. 

2.2.3 Potential Models and Predictions for New States 

Once the parameters of a specific potential model are fixed from the data on the known states of 
charmonium, it is natural to use the same approach for predicting the masses of yet unobserved res- 
onances corresponding to higher energy levels in the system. A large number of such predictions can 
be found in the literature spanning last three decades. Some of the results can be found in the papers 
[T5l IT8l IT9l |23| |32| [33| to name a few. The predictions for masses of the excited states considerably 
differ between models, which is not surprising in view of the nature of the approach. Moreover, given 
that the considered models do not have a controllable accuracy, it would be troublesome to assess what 
deviation of a prediction from the actual value of the mass should be regarded as successful. 

Nevertheless, such application of the models appears far from being an entirely empty exercise, and 
in a way provides some gross features of the expected spectrum of higher excitations. Namely all the 
'reasonable' models predict the same sequencing of levels: M{1D) < M{2P) < M[3S), which can 
possibly be traced to the general properties [20], and the specific masses being generally in the following 
ranges: M{1D) ^ 3.8-3.9 GeV, M(2P) ^ 3.9-4.0 GeV, and M{3S) > 4.0 GeV. The fine and hyperfine 
splittings of the levels are smaller than the uncertainty in the overall positions of the levels. All these 
states are above the DD threshold, so that most of them are expected to be broad due to decay into 
pairs of D mesons. The exception from this behavior can be found in the (2 ) and (2""'") 
resonances if they are below the DD* threshold at 3872 MeV. Indeed the unnatural spin-parity of these 
resonances forbids them to decay in the DD pairs, the only kinematically allowed states with open 
charm at such mass. The 2Pi and 3^5*0 states are expected to be well above the DD* threshold and for 
them such argument does not work. 

The issue of higher excitation of charmonium has recently gained a great attention due to observation 
of a whole 'zoo' of new charmonium-like states in experiment. In particular those, which seem to 
relatively well fit the expected pattern of the levels are the resonances Z(3930) and y(3940). The 
former state is found by Belle[35] in 77 production with the mass and width M{Z) = 3929 ± 5 ± 2 MeV 
and T{Z) = 29 ± 10 ± 2 MeV decaying mostly to DD, and which reasonably fits[19j the slot for 
the 2^P2 (Xc2) state of charmonium. The latter resonance, y(3940), observed by Belle|36] in the 
ujj/ip channel in the decays B —>■ uJ/ipK has the parameters M(Y) = 3943 ± 11 ± 13 MeV and 
T{Y) = 87 ±22 ±26 MeV. It appears to not decay into the pairs of pseudoscalar mesons DD and can be 
considered [21] as a candidate for the 2^Pi {x'ci) state of charmonium. If this interpretation proves to be 
correct an interesting spectroscopic question would arise related to the apparently inverted or small mass 
splitting between the 2^Pi and ^P2 states. Most recently the peak in the invariant mass of the system 
uiJ/%1) possibly consistent with y(3940) was observed by BaBar[37] in the decays ojJ/iIjK^ and 

B^ uoJ/ipKs- The values for the mass and width of the observed peak are somewhat off compared to 
the initial observation: M{Y) = 3914l| j ± 1.9 MeV and T{Y) = 33^1^ ± 5 MeV. Thus the status of this 
resonance is still not clear. The suggested interpretations include an excited P wave quarkonium[6|, 
hybrid ccg state|38j, and a four-quark molecular state. 

Another recently found resonance X{3940) is observed )40j as recoiling against J/ip in e+e^ 
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J/i/j + X, which imphes that its C parity is positive. Furthermore it appears to decay into DD* but not 
into DD, so that it hkely has unnatural spin-parity. These properties invite an interpretation [T^l EH] 
of the resonance as a O""*" charmonium state which would then be ?7c(3>S'). A possible problem of such 
interpretation is that most of the models expect the 35" level in charmonium to be somewhat above 
4.0 GeV, so that if further study of X{3940) indeed identifies it as a O""*" resonance, this may bring some 
new interesting understanding of hadron dynamics. 

2.3 Spectral Methods 

The potential models of heavy quarkonia and of charmonium in particular, are intuitively appealing, 
versatile and very convenient for estimates of various characteristics of the heavy resonances, but they 
can not be entirely satisfactory due to their model-dependent relation to the underlying theory of QCD. 
More directly related to the first principles of QCD are the methods based on the spectral relations 
for correlators in QCD. Such approach can be illustrated in its most basic form by considering the 
correlation function of the type F{x) = {0\T {O^{x),O{0)}\0), where 0{x) is a local operator and |0) 
is the vacuum state in QCD. Of relevance to charmonium is the choice of the operator 0{x) where it 
contains a factor cTc (with some structure F) and thus produces states of charmonium. The correlation 
function can be written in terms of the spectral sum over the physical states |n) containing a cc pair: 



where -D„(x) is the propagator of the state n. The lowest mass states contributing to the sum are the 
one-particle states i.e. the charmonium resonances, while at higher mass the sum is also contributed 
by the continuum of states containing the charm - anticharm quark pair. By an appropriate choice of 
the operator 0{x) one projects out the states with particular quantum numbers, while a suitable choice 
of X allows to make the sum being dominated by the charmonium resonances of interest. The direct 
relation to 'the first principles' arises when the correlator F can be also calculated in the interesting 
range of x by methods of the underlying QCD theory, thereby relating the phenomenological properties 
of hadrons to the results of a QCD calculation. 

The two approaches to calculating the correlator F are the numerical calculations in lattice QCD 
and the short-distance QCD analytical treatment. The lattice approach in principle allows to evaluate 
the correlator at large Euclidean separation x where the spectral sum is given by only the lowest mass 
state, so that e.g. the mass of this state can be fully determined. On the other hand, the short-distance 
QCD methods are restricted to relatively small values of the interval x, so that the spectral sum still 
contains some contribution from higher states as well as the lowest one. For this reason the relations 
resulting from calculations of this type are known as the QCD sum rules. The usefulness of the sum rules 
for description of the lowest state in a given channel depends on the existence of an intermediate range 
of a parameter analogous to x, where both the theoretical uncertainty in a short-distance calculation 
and the phenomenological uncertainty of the contribution of the higher mass states to the spectral sum 
can be reasonably controlled. 

2.3.1 QCD Sum Rules 

It is due to the e~^e~ annihilation data that the most well studied channel with hidden charm is the 
vector one, i.e. corresponding to the charm - anticharm production by the electromagnetic current 
of the charmed quarks = (c7^c), and this is the channel for which the original QCD sum rules 
were developed[41J. The relevant correlator is then the vacuum polarization P{q^)., considered in the 
momentum, rather than the position space: 



F{x) = Y.\{n\Om'D.^{x) 



(14) 



n 



P{q') {-q^9,u + q,qu)=t j d'xe^'^^ (0 |T {j^(a;), j.(0)}| 0) 



(15) 



10 



The spectral sum then takes the form of the dispersion relation for P{q^), 

TT J s [s — q'^ — le) 

and the imaginary part ImP(s) is related to the contribution of the states with hidden charm Rc{s) to 
the measured cross section ratio R{s) = cr(e+e^ hadrons) / {ATra"^ /3s): 

ImP(.) = ^ . (17) 

At values of sufficiently below the (perturbative) threshold at 4m^ the vacuum polarization P{q^) 
is determined by the QCD dynamics at short distances and can be calculated by using the Operator 
Product Expansion (OPE) for the T product in Eq. flT^ . 

T{j{x),3m = Y.u^)om (18) 

d 

in terms of local operators 0^(0) with increasing dimension d, and q(x) being the coefficient functions 
calculable in QCD. The leading operator of lowest dimension in this expansion is the unit operator 
X, and the corresponding coefficient Co(x) includes all the QCD perturbation theory result for P{q^). 
The first nontrivial operator in the series is the next one with dimension d = A and is quadratic in the 
gluon field strength tensor G^^,, so that its contribution to the vacuum polarization is proportional to 
the gluon vacuum condensate |42[ HS], HI HI]. The relation between the theoretical expression for P{q^) 
and the phenomenological integral over the observed cross section can be studied as a function of q^ far 
below the threshold. Alternatively, one can compare the expressions for the derivatives of the vacuum 
polarization with respect to at = 0, for which the 'phenomenological' side of the sum rules is given 
by the moments of the ratio Re'. 

■^n= f ^ds, n=l,2,... , (19) 

J s 

J\A.n = (127r^/n!) {d^P{z) / dz'^)\z=q2=Q. Theoretically, the Taylor expansion of the vacuum polarization 
including the leading nonperturbative term can be written as 

where C„ and P„ are dimensionless coefficients that are calculated as series in powers of a^- The 
coefficients C„ are known [45J at arbitrary n up to a^, and Ci has been recently evaluated P6| HT] to 
order a^. The coefficients P„ are known in the lowest [121, H] and the next to lowest [48j orders in a^. 
The contribution of a narrow resonance to Rc can be approximated by a 5 function: 

R,[s) = —5{s-M^)T,,M, (21) 



a 



where M is the mass of the resonance and Fee is the width of its decay into e^e~ . The lowest charmonium 
state contributing to the dispersion integral in Eq. (fT6|) is the J / ip resonance, so that the n— th moment 
can be written as 



Mn 



ree(J/V') ^ ree(V'') 



+ / ^ds. (22) 



where the latter integral runs over cm. energies above the resonance. One can readily see that 
the relative weight of the lowest resonance grows with the number of the moment n. (E.g. the J/V' 
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contribution essentially dominates the moments of charm cross section already at = 3 -i- 4.) However 
the correction terms, both the perturbative and nonperturbative, in the theoretical calculation of the 
moments also grow with n, which is a reflection of the fact that higher derivatives are sensitive to the 
threshold singularity in the correlator where the essential distances are no longer short. In particular 
the parameter for the perturbative expansion for the moments is effectively as \fn (corresponding to 
the agjv behavior of effects of the Coulomb-like interaction near the threshold), and for the leading 
nonperturbative term due to the gluon condensate behaves as r? jm^. It was found nevertheless |4T| 
H] that the perturbative one-loop expression for the moments can be relied on up to n f» 4, and somewhat 
higher moments (up to n = 7-r-8) then can be used in order to determine |42 1 H3] the value of the gluon 
vacuum condensate (G^). Once other relevant QCD parameters (a^, the charmed quark mass 
were determined from the data and the sum rules for the vector channel, they could be used for other 
channels as well. In particular, the sum rules for the spectral density of the pseudoscalar operator (0750) 
correctly predicted [42j| the mass of the lowest 0""'' charmonium resonance, the ?7c. 

The best theoretically known is the first moment of Rc{s), due to the accuracy (a^) of the available 
perturbative calculation and due to very small nonperturbative term proportional to (G^). On the 
experimental side this moment is rather sensitive to the details of the charm production cross section in 
e~^e~ annihilation at around 4.0 GeV where the open charm production sets in, and the energy behavior 
of the cross section is quite complicated. Recent data in this region [49t |50] have allowed to put to use 
the attained theoretical accuracy in A^i. In recent analyses the sum rule for the first charm moment 
was used for precision determination of the short-distance charmed quark mass parameter, Wicirn^: 
1295 ± 15 MeV[i7] and 1286 ± 13 MeV[5T]. 

2.3.2 Lattice Methods and Limitations of the Spectral Approach 

The QCD sum rule approach, based on the OPE and analytical calculations is certainly limited by the 
applicability of the perturbative expansion in as- This requires a careful choice of parameters, e.g. the 
number of the moment n, to ensure such applicability and still get a phenomenologically useful relation. 
It is widely believed that a ticket to calculations beyond the perturbation theory in ag that are not 
limited to short-distance QCD is the numerical lattice approach. In particular, a natural application 
of this approach would be a calculation of correlators at large separation between the points in the 
Euclidean space, so that the contribution of the lowest states in each channel could be completely 
determined. However practical implementations of lattice calculations to heavy quarkonium including 
charmonium still run into difficulties of their own. A detailed discussion of recent developments of 
the lattice methods and of the associated difficulties can be found in the review[2]. As of this writing 
the accuracy of the lattice results even for the masses of the lowest states of charmonium in each J^*-^ 
channel leaves an ample room for further improvement. 

The spectral approach in general, based on analytical or numerical calculations, has 'built in' a 
certain deficiency with regards to excited states. Namely the spectral sum for a correlator at a Euclidean 
separation necessarily receives the largest contribution from the lowest physical state. By increasing 
the separation one can at best enhance the sensitivity to the properties of the lowest state and thereby 
evaluate those properties. There is however no simple way of 'focusing' spectral relations on excited 
states, so that a study of those states by spectral methods runs into additional difficulties. 

Furthermore, even at the level of sum rules, or of a study of the lowest charmonium states in each 
channel by lattice methods, the spectral approach eventually runs into limits of its accuracy arising 
due to annihilation of charmonium into light hadrons. This limitation can also be viewed as a version 
of the problem of excited states: the lowest state of charmonium in a channel with given quantum 
numbers J^*" is certainly not the lowest hadronic state in that channel. In other words, at some level 
an operator (cFc) produces states containing only light quarks and gluons and no charmed quarks. The 
old phenomenological rule, according to which a mixing of hidden quark flavors, although not forbidden 
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a) 



Figure 2: Representative graphs for OZI rule violation in the vacuum polarization by the elec- 
tromagnetic current jc of charmed quarks (a) and for the interference between the charmed 
quark current and that of the light quarks (6). The wavy lines show the gluons. 



by any conservation laws, is dynamically suppressed, is traditionally referred to as the Okubo[52]- 
Zweig[53]-Iizuka[M] (OZI) rule, so that the discussed effect is precisely the violation of this rule. The 
significance of this effect generally depends on the channel considered, and as also will be mentioned in 
the Section 3.2.3, it can be essential in the properties of the rjc resonance. In the considered above case of 
the vector channel this effect is likely to be quite small, fully in line with the known strong suppression 
of OZI violation in vector mesons. Indeed, the partial width of the J/ip resonance associated with 
strong annihilation of the hidden charm into light hadrons is only about 70 keV, which can be viewed 
as the imaginary part of the shift of the mass of the resonance due to the OZI violation. Generally the 
real part of the shift is of the same order as the imaginary, which for the J/ijj meson would be much 
smaller than any current accuracy of a theoretical calculation of its mass. 

It can however be noted that the already achieved accuracy in the first moment of Rc{s) is actually 
on the verge of being sensitive to the effects of the OZI rule violation. Indeed, the separation of the 
hidden charm from light degrees of freedom in the electromagnetic production is broken in order due 
to the mechanisms illustrated in Fig 121 The theoretical calculation of the first moment Aii happens to 
be insensitive |55[ [55] to the contribution of the three-gluon intermediate states shown in Fig 12^, while 
the cross term between the currents of the light and charmed quarks, shown in FigJ2)D does not enter the 
correlator in Eq. f|T5|) by definition. The absence of the contribution from the graphs of the type shown 
in Figl2^ in A^i (as well as in Ai2 and Ais) can be assured by the following reasoning. Consider this 
mechanism at small g^, so that <C m^. In this region the charmed quark loop can be contracted into 
a point thus generating an effective point-like Lagrangian for coupling of the current to three gluons. 
Similarly to the well-known Heisenberg-Euler Lagrangian for the photons, the current conservation and 
the QCD gauge invariance ensure that the mass of the fermion in the loop enters as l/m^. The graph 
of Figl2^ contains two such heavy quark loops, so that the low-energy expansion of the contribution of 
the discussed mechanism to the correlator of the currents starts as {q^ /fnlY ^"^{(f / ^V) ■> therefore 
no contribution to the first three moments M. results from this mechanism. 

The cross-talk between the electromagnetic currents of the light and charmed quarks illustrated in 
Figl2]D in fact relates to the general problem of phenomenological separation of the hidden-charm and 
light-quark contribution both to the measured cross section and to the sources of the observed final 
states. The former separation of the hidden-charm part from the total cross section is usually done 
by subtracting from the data the smooth background associated with the light states, while the latter 
separation of the sources corresponds to evaluating how much of light states are produced by the current 
of charmed quarks and how much of the hidden charm is produced by the electromagnetic current of the 
light quarks. The cross-talk effect, however is suppressed by the flavor SUidt) symmetry. Indeed, the 
electromagnetic current of the light quarks is a pure component of SU{?i) fi octet, so that in the limit of 
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exact symmetry this current does not produce the heavy quarks, which are SU (3) fi singlets. This effect 
is further suppressed in the first moment Aii (in this case only in the first moment). The reasoning 
is similar to the one regarding Figl2^. In this case there is only one heavy quark loop proportional 
to l/m'^. Thus the low-energy expansion of the total sum of the graphs of the type shown in Figl2)D 
necessarily starts with {q^ /mIY giving no contribution to M.i. 

It can be also mentioned that in the vector channel the mixing of light hadrons and the hidden 
charm starts in the order due to the C parity. In the C-even channels such mixing starts in the 
order (except for the J = 1 channels, where it is forbidden by the Landau- Yang theorem [57[ 1581 H]). 
so that the OZI violating effects should be larger in the spectroscopy of these channels. Furthermore, 
the J = channels O""*" and O"^"*" generally suffer from large quarks-glue mixing effects, related to direct 
instantons[59], which may hold the clues to understanding some properties of the r/c and Xo resonances. 

3 Charmonium Annihilation 

The same OZI violating effects that are complicating an improvement in the precision of calculation 
of the masses of charmonium resonances are also the sole reason for the eventual disappearance of 
the hidden charm, i.e. for annihilation of the cc quark pair into light states. The two interactions 
contributing to the annihilation decays are the strong interaction responsible for most processes of this 
type, and the electromagnetic interaction which is relevant to annihilation decays of the vector 1 
states and the 27 decay widths of the C even states with J 7^ 1. We first consider in this section the 
electromagnetic processes. 

3.1 Electromagnetic Annihilation 

3.1.1 Annihilation of^Si States through Virtual Photon 

The ^5*1 states have quantum numbers of a virtual photon, J^^ = 1 and can annihilate into lepton 
pairs or light hadrons through one photon. This is also the process, which when reversed, gives rise to the 
formation of the ^Si states as resonances in e~^e~ annihilation. The rate of the decay can be estimated 
in the extreme-nonrelativistic picture, where the system is described by the wave function ip{r) for the 
quark- ant iquark pair and depending on their relative position f = fc — r^. The annihilation takes place 
at the characteristic distances of order 1 / which are to be viewed as r — > for a nonrelativistic pair, 
so that the annihilation amplitude is proportional to the wave function at the origin. For an n-th S-wave 
state the wave function can be written entirely in terms of its radial part Rns'- i^nsi^ = Rns/V^, and 
the expression for the rate of annihilation into e^e~ takes the form 

rUn'S,) = |i?„,(0)p (1 - ^) , (23) 

where Cc = 2/3 is the electric charge of the charmed quark in units of the fundamental charge, M is the 
mass of charmonium, and finally, the term with in the parenthesis gives the first QCD correction. 
The coupling constants and a should be taken at the scale rric (or M) as appropriate for a process 
proceeding at distances ~ l/rric. The correction (sometimes forgotten) due to the running of the QED 
coupling is reasonably tractable and results in approximately 7% enhancement of the rate, while the first 
QCD correction is quite discouraging. Indeed its numerical value, — l.Ta^, amounts to —(0.35 -7- 0.5) 
for as{mc) in the realistic range 0.2 -^ 0.3, which tells us that higher QCD corrections can be quite 
essential [60]. 

There are however other uncertainties of a conceptual nature associated with the formula (123|) . One 
source of such uncertainties arises from the use of the nonrelativistic (essentially static) approximation. 
One of the indications of this approximation is the parameter M: "the mass of charmonium" entering 
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Eq. fl23l) . It would be impossible to specify at this level of approximation whether M should be set at 
twice the quark mass 2mc, or the mass of the specific decaying state, or some combination of those, 
since the difference between these values for M is formally of order f^/c^. Clearly, the first relativistic 
correction can amount to tens percent. It should be noted that the specific form of this correction does 
not come from purely kinematical effects in the annihilation of moving, as opposed to static, quarks, 
but is also sensitive to the interaction, since the average kinetic energy is of the order of the average 
interaction potential as follows from the virial theorem. Another fundamental uncertainty of Eq. (!23ll 
is related to the assumption that the quarkonium can be described by a quark-antiquark two-body 
wave function, which assumption is formally not valid in nonperturbative QCD, as was discussed in 
Sec. 2. 2.1. The nonperturbative correction to the width Fee of the ^5*1 states has been calculated |13] 
in the limit of very heavy quarkonium, where the correction is parametrically small. An extrapolation 
down to charmonium would produce an unreasonably large result [60j, neither it would be justified. 
Thus the issue of the applicability of a description of charmonium as a two-body system arises again 
as previously in the general discussion of potential models. If in spite of this issue potential models are 
used in conjunction with Eq. (123!) and similar formulas for other annihilation rates that will be discussed 
in this section, the results are in a qualitative agreement with the observed pattern of the decay rates, 
although no controllable accuracy can be assigned to such results. 

The electromagnetic annihilation contributes a sizable fraction, about a quarter, of the total decay 
width of the J /if) resonance. The latest precise data|6T] have moved the world average [27] for the e^e~ 
branching ratio to BeeiJ/'ip) = (5.94±0.06)%. The decay into /i"*"yU~ goes with essentially the same rate, 
and the rate of electromagnetic annihilation into light hadrons is given by the e~^e~ rate scaled by the 
ratio R measured in e'^e~ annihilation just below the J/ip resonance |49j: R{3.0 GeV) = 2.21±0.05±0.11, 
thus giving the total branching fraction for the decays of J/ip through a virtual photon as 

B{J/^ ^ 7* X) = (2 + i?) Bee{J/tp) = (25.0 ± 0.8)% . (24) 

For the ip' resonance the contribution of the electromagnetic annihilation is by far less prominent 
(partly due to larger total decay width), given that Bedi'') = (7.35 ± 0.18) x 10^^ and using[l9] 
i?(3.7GeV) = 2.23 ± 0.08 ± 0.08, one can estimate 

B{ip' ^ 7* ^ X) = (2 + 0.39 + R) i3ee(V^') = (3.39 ± 0.11)% , (25) 

where the term 0.39 in the parenthesis accounts for the decay ip' r+r^. The theoretical ratio of the 
rate of this decay to Teeii'') is in agreement with the direct measurement [27], but has about ten times 
smaller error. In terms of potential models the ratio of the e~^e~ decay rates for ip' and J/ip can serve 
according to Eq. (l2^ as a measure of the ratio of the wave functions at the origin, and this estimate 
was used above in connection with the relation in Eq. ffTOl) . 

The e'^e~ decay rates of the vector resonances appear directly in the spectral approach (Eq. fl^ ). so 
that the relativistic effects and perturbative and nonperturbative QCD terms are all included inasmuch 
as these effects are taken into account in calculation of the correlation function. However, as mentioned 
previously, the spectral relations are mostly sensitive to the contribution of the lowest resonance, i.e. 
the J/ip in charmonium. It is known since long ago [HI H] that the QCD sum rules are in an excellent 
agreement with the observed value of Fee(J/V') of about 5keV. 

3.1.2 Annihilation of mixed ^Di — ■^S'l States through Virtual Photon 

The quantum numbers of a ^Di state are J^*^ = 1 and allow such state to decay through one virtual 
photon. However in a pure D wave state the wave function at the origin is vanishing, so that in the 
leading nonrelativistic approximation the annihilation amplitude is zero. A non-vanishing contribution 
to this amplitude arises in the order v"^/ due to two mechanisms. The first one is the ^Di Si mixing 
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due to the tensor force, while the second arises from the expansion of the D wave coordinate wave 
function at the annihilation distances 1 / rric and is proportional to the second derivative of the radial 
function at the origin: i?"^(0)/m^, which is also of order /w? ~ f^/c^ in comparison with the S'-wave 
wave function at the origin. The latter mechanism alone would result in the expression for the width [1] 

Tee{n'D,) = '^^^\Knm\ (26) 

where no short-distance QCD correction is included. The direct annihilation amplitude and the one 
due to the mixing fully interfere with each other and are of the same order in v'^/c?. Thus it would 
generally be unjustified to consider one mechanism but not the other. 

Phenomenologically, the resonance V'lSTTO) is considered to be such dominantly l^Di state with an 
admixture of ^5*1 wave function. The simplest model for the latter admixture uses the proximity in mass 
of the ip' resonance and considers only the two state '?/'(3770) — ip' mixture. Then, using the data[27] for 

ee 

[^(3770)] = 0.242t o;o24keV and ignoring the direct annihilation amplitude one can estimate [62l [63] 
the mixing angle 9 as being about 0.2. Such estimate certainly agrees with the expectation for the size 
of the f^/c^ effects in charmonium, although the particular numerical value should likely be taken with 
certain reservations given the very simplistic nature of the model. 



3.1.3 Two-photon Annihilation of C -even States 



The cc quark pair in C even states with J 7^ 1 can annihilate into two photons [64l H]. For the u^Sq 
states the amplitude is proportional, in a potential model approach, to the wave function at the origin, 
-Rns(O), SO that it makes sense[65j to consider the ratio of the ^Sq 27 and J/ip — >• e+e~ decay rates 
where the value of the wave function at the origin cancels. Including also the first short- distance QCD 
correction [66j for the 27 decay, which can be traced back to the positronium result [67], one can write 
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(27) 



Experimentally the 27 decay rate is measured for the t^c, albeit with a large uncertainty [2 7], so that for 
the 15 charmonium the ratio of the rates in Eq. fl27|) is experimentally 1.3 ± 0.4. Although this value 
can be considered as being in agreement with the theoretical expectation, measurements with smaller 
experimental uncertainty are clearly needed for a more meaningful comparison. It would also be of a 
great interest to test Eq. fl7r|) for the 25 charmonium, i.e. for the corresponding decay rates of the tp' 
and T]'^ resonances. So far the process 27 — >• 77^ has been seen [681 EH] as a source of rj'^ in 77 collisions, 
however quantitative data on decay widths of the 77^ are still in flux. 

It can be mentioned that the rate of the decay 77^ 27 can be analysed by considering the QCD 
sum rules for the moments of the 77 scattering cross section due to the electromagnetic current of 
the charmed quarks [701 H]. The rj^ resonance is the lowest state contributing to this cross section for 
perpendicular linear polarizations of the photons. The result of such estimate, not including QCD 
corrections and nonperturbative terms is[4J r(?7c ~^ 27) = (6.1 ^ 6.5) keV, which is also in agreement 
with the existing data[27]. 

The amplitudes of the two-photon annihilation of the P-wave states, ^Pq and ^P2 are proportional 
to the first derivative of the radial wave function at the origin: R'p{0)/mc. The specific expressions[ 
with the first short-distance QCD correction [TTj are given by 
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Both of these rates are of order jc? in comparison with the similar rate for the ^5*0 state. Thus one 
should expect a certain suppression of these decays of the Xcj resonances as compared to the rate V{r]f. 
77), although an absolute prediction of the rates would be quite model dependent due to uncertainty 
in the value of the wave function (as well as due to other uncertain factors discussed in connection 
with Eq. (l23l) ). The experimental values [27] for these decay rates are: r(xco 11) = 2.9 ± 0.4keV 
and r(xc2 11) = 0.534 ± 0.050 keV, which can indeed be considered as somewhat suppressed in 
comparison with T{ric —>■ 77). As usual, some of the theoretical uncertainty goes away if one considers 
an appropriate ratio, in this case the ratio of the two decay rates in Eq. (!28l) : 



r(3p, ^ 77) _ 4 
r(3Po ^ 77) 15 



37r V 3 



±(1-5,51^1 , (29) 



The experimental value of the ratio calculated from the world averages|27j for the decay rates, 0.185 ± 
0.025 is indeed smaller than the uncorrected value 4/15 = 0.267 in Eq. (l29|) . which agrees with the 
negative value of the QCD correction. However, the latest data with the best precision in a single 
experiment [?] give 0.235 ± 0.042 ± 0.005 ± 0.030. In either case it would still be troublesome to draw a 
more quantitative conclusion due to the large coefficient of the first radiative term. 

The same considerations regarding the annihilation to 77 can be applied to the recently observed 
state Z(3930) interpreted as the radial excitation Xc2- The observation of this resonance |35j is in fact 
due to its discussed coupling to 27. However no quantitative data on the electromagnetic decay width 
are available so far. 

The subject of the two-photon annihilation may also become of relevance for the yet unobserved 
1^D2 (2""'") state of charmonium. If the mass of this resonance is below the DD* threshold it should 
be quite narrow since its unnatural spin-parity forbids decay into DD pairs. No mixing due to the 
tensor forces is possible for this state, therefore its two-photon annihilation would provide an access to 
the amplitude of direct annihilation from D-wave. A similar amplitude enters the previously discussed 
annihilation of a ^Di state into e~^e~, where however it gets tangled with the S — D mixing effects. 
Thus a measurement of the decay 1^D2 77 can be helpful in separating the direct annihilation from 
mixing in the properties of '?/'(3770). The rate of the decay is given[4j as 

2^ 3p^n'2 

T{n'D2 - 77) = — ^ \Knm' (30) 
in terms of the second derivative at the origin of the same radial function as in Eq. 



3.1.4 J/tp 777 

The decays ^5*1 —>■ 37 have very small rates proportional to a^. However a measurement of such decay 
for the J/ip resonance does not appear to be unrealistic. The rate of the decay in the potential approach 
is given by 

^ / T / , X 16 (vr^ — 9) , , mo / as\ , , 

TiJ/^ - 37) = ^ 3^ J - \Rism' (1 - 12.6^) , (31) 

where the lowest-order result is an adaptation of the orthopositronium decay formula [73] and the first 
QCD correction is also an adaptation [65] of the numerical result [74] for the one-loop QED correction to 
the orthopositronium decay rate. As usual, model-dependence can be reduced by considering the ratio 
of this decay rate to Tee (Eq. (l23|) ): 



3-,) _ 64 (.^ - 9) ^ . _ ^ gO,-. ^ ^ 3 ^ . _ ^ 



Tee{J/ip) 2437r V ■ vr / ■ V ' ir 

The QCD radiative correction in this estimate is large, and it is not clear what numerical value it 
should be assigned. The uncorrected number puts the three-photon decay rate in the ballpark of 3 eV 
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corresponding to the branching fraction B{J/il) — > 87) ~ 3 x 10~^, which can be compared with the 
current upper hmit (at 90% CL): 5.5 x 10"^. 

The three-photon decay, inspite of its small rate, presents a very interesting object for a study 
of dynamics of charmonium. Indeed, the ratio of the rates in Eg. (132!) is sensitive to only the QCD 
corrections, so that a measurement of this ratio can possibly shed some light on understanding of the 
behavior of the QCD radiative effects in a situation where the coefficients in the loop expansion are 
large. Since such behavior of the expansion coefficients is typical for various effects in heavy quarkonium 
a better understanding of the QCD expansion can be of help in considering those other effects as well. 
Furthermore, it would be of great interest if one-photon energy spectrum in the 87 decay could be 
studied, since the photons provide a very clean "CT scan" of the internal structure of J/ip. Namely, 
the photon spectrum at energy uj provides sensitivity to distances ~ 1/ ^JirtcOJ, so that dynamics of the 
charmed quarks can be probed by this spectrum at distances ranging from the short ones, ~ I/wt-c, up 
to the typical charmonium size. 

3.2 Strong Annihilation into Light Hadrons 

Charmonium decay into light hadrons through the strong interaction is viewed in QCD as a two-stage 
factorized process. First the cc pair annihilates into gluons at 'short' distances of order l/rric and 
then the gluons fragment into specific hadronic final states. The total inclusive probability of the 
latter fragmentation is considered to be equal to one, so that the total decay rate is determined by 
the short- distance annihilation to on-shell gluons. Clearly, this approach based on a perfect 'gluon - 
hadron duality', involves an uncertainty of its own in addition to the previously mentioned uncertainties 
involved in calculation of charmonium annihilation. Inspite of this reservation, such treatment first 
considered in Ref. [75], is extremely successful in explaining, at least semi-quantitatively, the pattern 
of OZI violating narrow widths of the charmonium resonances below the open charm threshold, in 
particular in understanding the very narrow width of the J / ip resonance, which has so profoundly awed 
particle physicists in November 1974. 

3.2.1 Three- gluon Annihilation of ^ Si Charmonium. 

The minimal number of gluons into which a ^Si state of a heavy quark pair can annihilate is three, since 
the process through one virtual gluon is forbidden by color and a two-gluon final state is excluded by 
the negative C parity of the initial state. The decay rate in the lowest order in QCD [75] can be found 
by decorating the orthopositronium decay formula [73] with the appropriate color factor, while the first 
QCD radiative correction is known only numerically [76], and being expressed in terms of as normalized 
at the scale rric in the MS scheme reads as 

Tin's, ^3g^ light hadrons) = ^ ^ |i?.,(0)r fl - 3.7^) . (33) 

The ratio of this rate to Fee is then sensitive, in this approach, to only the coupling constants: 
T{n'Si ^3g ^ light hadrons) 5 tt^ - 9 a^(mc) A , „ «s\ 

ree(^ Oi) 18 TT a'^ \ TT / 

One is naturally tempted to use this formula for an estimate of the QCD coupling as- In order to 
make such estimate one has to evaluate the branching fraction for the direct strong annihilation of 
J/ip, which amounts to subtraction from the total sum the electromagnetic annihilation contribution, 
estimated in Eq. (!24|) . the small contribution of the decay J/ip ^ 7??c with the branching fraction [27] 
(1.3 ±0.4)%, and the contribution of the direct photon emission, originating in the process J/ip — > 7(7(7, 
which will be discussed in some detail few lines below. The branching fraction for the latter process is 
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somewhat uncertain experimentally due to a large background of secondary photons in the low-energy 
part of the photon spectrum. Namely, the direct photon emission has been observed and reliably 
measured |77j at x > 0.6, where x = 2uj^/Mj/^ is the ratio of the photon energy to the maximally 
allowed by kinematics. The integral over the observed spectrum in this region corresponds to the 
branching fraction B^J/ip — > igg)\x>os = (4.1 ± 0.8)%. The experimental shape of the spectrum 
reasonably suggests that the observed part makes about one half of the total decay rate with another 
half being hidden under the background at x < 0.6. Thus one can rather conservatively estimate the 
total branching fraction as B{J/tp — > jgg) ~ (8 ± 3)%. As a result the fraction of the direct strong 
annihilation of J/ijj can be estimated as l3{J/ijj 3g) ~ (66 ± 3)%, and the ratio in Eq. llMl) a^ 

5<#^ = 11.1 ±0.5, (35) 

thus providing one with the estimate as{mc) ~ 0.19. It would however be difficult to assign a reliable 
error bar to this number, given the reservations mentioned previously regarding the assumptions and 
approximations made in connection with the described calculations of the annihilation rates. 

A marginal suitability of the direct hadronic decay rate of charmonium for precision determination 
of QCD parameters can also be illustrated by comparing the data on the annihilation decays of J/tp and 
The proportionality of the annihilation rates to \Rns{0)\'^ implies a similarity between the decays 
of J/ip and ip'- Namely, the ratios between these two resonances of their similar decay rates should be 
the all equal to each other. In particular the ratia of the branching fractions for similar decays should 
all be equal to Beeiip')/B{J/'ilj) = (12.4 ±0.3)% (the so-called "12% rule"). A recent CLEO-c dedicated 
study[78] of the decays of ip' ending up in J/ip in the final state allowes to separate the branching 
fraction for the direct decays of ip' into light hadrons, which includes the electromagnetic decays of 
this type and the 'ygg decays with direct photon: B{ip' — *■ lighthadrons) = (16.9 ± 2.6)%. The same 
branching fraction for J/tp is listed in the Tables [27]: B{J/ip —>■ lighthadrons) = (87.7 ± 0.5)%. The 
ratio of these numbers gives B{ip' — lighthadrons)/;B(J/'?/' —>■ lighthadrons) = (19.3 ± 3.0)%, which 
substantially differs from the "12% rule". One would also arrive at a similar contradiction with the 
simple picture of charmonium annihilation if first subtracted from B{tp' — > lighthadrons) the small 
contribution of the electromagnetic processes ip' —>■ 'j* —>■ lighthadrons and the (estimated) fraction of 
ip' — > 'jgg, leaving B{ip' — > 3^^) ~ (14.3 ± 2.9)%, which gives the estimate for the ratio of the ggg and 
e'^e~ annihilation rates of ip' equal to 19.5 ± 4.0, i.e. different from the same ratio for J/ip in Eq. (l35p . 

That there is more dynamics to the hadronic annihilation of the ^Si charmonium than the simple 
factorized approach suggests, is indicated by that the 12% rule is also quite conspicuously broken in 
exclusive decay modes, the most famous in this respect being the decay into pvr which is strongly 
suppressed for ip', well below the "12% rule". A number of other exclusive modes with a (seemingly 
random) deviation from this rule have been measured experiment ally [271 EH]. Although deviations from 
the simple picture are expected on general grounds, the specific mechanisms for such deviations are 
currently unknown even at the level of inclusive decay rates. We will discuss some considerations on 
this subject at the end of the current section on charmonium annihilation. 

3.2.2 The Decay into 'jgg 

The mixed electromagnetic and strong annihilation into a photon and two gluons, briefiy discussed 
above in connection with the estimate in Eq. (l35p . is actually quite interesting on its own[SIJ]. The ratio 
of the total decay rate to that of the decay into e~^e~ in the simple annihilation picture depends only 
on the QED and QCD coupling constants: 

T{n^Si^-fgg) 8 tt^ - 9 ^^(^J / « 



Teein'^Si) 9 TT a 



1.3— , (36) 

TT 



^This estimate also agrees with the result of a recent similar evaluation[6j done in a slightly different way. 
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where the QCD radiative correction is from the resuh of numerical calculation of Ref. [76]. Using the 
known value of BeeiJ/'ip) and ^^(mc) ~ 0.19 as estimated from Eq. flS^ . one can estimate B{J/ip —>■ 
^gg) ~ 6.7%, which is in a reasonable agreement with the available data[77j. 

Besides the total rate of this radiative hadronic decay, of a great interest is the spectrum of the direct 
photons[4], since this spectrum provides some insight into how the gluon - hadron duality sets in. Indeed 
the invariant mass squared, g^, of the hadronic final state into which the two gluons fragment is related to 
the photon energy u; as = Mjy^ (1 — x), where x = 2u;^/Mj/^. The 'parton' spectrum, corresponding 
to on-shell gluons extends all the way to the kinematical boundary at x = 1, corresponding to = 0. 
In reality, one certainly does not expect a gluon-hadron duality at low g^, so that the spectrum of direct 
photons at high energy should be suppressed by hadronic effects, as is observed in experiment [77]. At 
larger g^, corresponding to lower x, the actual spectrum should approach the parton curve due to the 
onset of the gluon-hadron duality. It should be noted however that, as previously mentioned, the direct 
photon spectrum at lower x (in practice at x ~ 0.5 for charmonium) goes under an overwhelming 
background of secondary photons which are present in hadronic final states. Thus in the decays of 
charmonium it is unlikely that the fragmentation can be studied at g^ beyond 4 or 5 GeV^, and a much 
better testing ground for such study is provided by the decays of the bottomonium T resonances [81]. 



3.2.3 Two-gluon Annihilation of C -even States 

The C-even states of quarkonium with J 7^ 1 can annihilate into two gluons, much in the same way 
as they decay into two photons. In fact in the lowest order the gg and 77 decay rates are simply 
related [611 H] as F^g/r^^ = 2ag/(9e^Q;^) = (9/8) (tts/a)^. The coefficients of the first QCD correction 
are, naturally, different for individual states [66l [71] : 



99) 9 



r(^^o - 77) 



a 



1 + 8.2- 



(37) 



r(3Po ^ gg) 9 



rppn 



77) 



a 



n 2 



1 + 9.3— 

TT 



r(3p2 ^ gg) 9 



77) 



a 



n 2 



1 + 3.1 — 

TT 
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The correction for the decay rate of the state is presently unknown. One can see that the coefficients 
of the first QCD correction for ^5*0 and ^Pq states are unusually large, which makes difficult a reliable 
quantitative comparison with the data. The uncorrected formulas with ~ 0.2 would give for the ratio 
of the width the value of about 850, which is not even close to the observed values[2^: V gg{rjc) /T ^^(jj^ = 
(3.57+1.15) X 103, r,,(xco)/r^7(Xco) = (3.62 + 0.43) x 10^, and Tggixc2) /T^^{xc2) = (3.09 + 0.23) x lO^. 
The experimental data can be somewhat understood using a significantly larger QCD coupling as{mc) ~ 
0.3 than is inferred from the three-gluon decays of J/ip. Such larger value of as{mc) better complies 
with the determination [22] of the QCD coupling from r decays and e+e~ annihilation data, but then it 
entirely misses the observed rate of the three-gluon annihilation. It is not clear at present how to avoid 
this crunch. One could possibly argue for the rjc and Xco resonances that their hadronic decay rates 
are enhanced by the instanton-type quarkonium glue mixing in the ^ and O"*"^ channels ||59j. However 
such argument would still leave the observed 2g/2'j ratio for the Xc2 unexplained. 



3.2.4 Strong Annihilation of ^Pi and ^Pi States 

The so far discussed amplitudes of the decays into on-shell gluons are not sensitive to the infrared 
behavior of quarks and gluons, at least at the considered here level of perturbative calculation at the 
leading and the one-loop level. In other words, these amplitudes are all expressed in terms of the wave 
function of the quark pair at the typical annihilation distances of order l/rric. In the two-gluon decays 
this behavior is guaranteed by the fixed and large energy of each of the gluons in the cm. system. 
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while for the three-gluon decays of the ^5*1 states the kinematical region, where one of the gluons is 
soft, is not enhanced due to considering the heavy quarks being essentially static and thus not radiating 
soft gluons. Clearly, such "good" infrared behavior generally becomes invalid in higher orders of the 
perturbative expansion, as well as beyond the leading nonrelativistic approximation. For the ^Pi and 
^Pi states the infrared behavior of the amplitudes describing their annihilation shows up already in the 
leading approximation [83]. The reason for the infrared sensitivity of these amplitudes is that a P wave 
necessarily involves a motion of the quarks, so that a P-wave state can go into an S'-wave by an El 
emission of a soft gluon. The spin of the quarks is not involved in the emission. Therefore the colorless 
^Pi state goes into a color-octet ^5*0 state of the quark pair, while the ^Pi state goes into a color-octet 
^Si. The intermediate ^5*0 state decays into two gluons, while the colored ^Si can decay either into 
a light quark-antiquark pair, qq through a virtual gluon, or, generally, into a gluon pair. It turns out 
however [83], H] that for a static colored ^5*1 pair of heavy quarks the amplitude of annihilation into two 
gluons is identically zero in the leading order, thus the resulting decay processes are ^Pi — >• ggg and 
^Pi —>■ gqq. In calculating the total probability one has to integrate over the energy uj of the soft gluon, 
which integration, as is standard for a soft gauge quantum emission, is infrared- divergent as / duj/uo. 
With a logarithmic accuracy the upper limit in this integral is set by the quark mass, and the lower 
limit is a typical QCD mass scale. 

The leading-logarithm expressions for the hadronic widths of the ^Pi and ^Pi states are given as [551 

TCP, ^3g^ hadrons) = ^ ^ |i?i>(0)P In , (39) 

9 Aqcd 

and 

Ti^Pi hadrons) = V Ti^Pi^gqq^ hadrons) = ^ |i?UO)|Mn . (40) 

One can see that the numerical coefficients in these formulas correspond to the hadronic width of the 
^Pi being larger than that of the ^Pi by the factor Svr/G ~ 2.6. The experimental value of the hadronic 
width of the l^Pi charmonium resonance, Xci, is approximately 0.6 MeV, so that based on these formulas 
one would expect the hadronic width of the l^Pi resonance, he, to be about 1.5 MeV. The latter value 
does not contradict the CLEO-c data[30j, but it might be in a contradiction with the Fermilab E835 
results [29] claiming the width of the observed he resonance to be F < 1 MeV. Apriori one would expect 
the applicability of these logarithmic estimates of the annihilation rates to be marginal at best. 



3.3 Non-perturbative Effects in Hadronic Annihilation 

The logarithmic formulas become formally applicable in the limit of very heavy quarks, where the 
resulting logarithm can be considered as a large parameter, and one can develop a formalism|84] of 
a logarithmic theory of the processes with soft-gluon transitions between color-octet and color-singlet 
heavy quarkonium states, or consider purely phenomenologically the mixing between pure colorless 
cc states and states where a color-octet quark pair is accompanied by a soft gluon field. An exten- 
sive overview of this approach, called in the literature Non-Relativistic QCD (NRQCD), and a list of 
references can be found in the review [5]. 

The effects of soft gluons can be addressed not only in the decays of the ^Pi and ^Pi states, where 
such effects are dominant, but also for 'well behaved' processes, where these effects give rise to small 
or moderate corrections. In particular the corrections arising from nonperturbative soft gluon field 
in the annihilation of the ^5*1 states may be responsible for the previously discussed deviations from 
predictions based on the simple perturbative picture, such as the "12% rule" . A consistent analysis of 
the leading nonperturbative corrections to the annihilation rates is however possible only in the limit 
of very heavy quarkonium, for which such corrections are expressed [85] in terms of the gluon vacuum 
condensate. 
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Figure 3: The mechanism describing the leading nonperturbative correction to the annihi- 
lation rate of a heavy ^Si quarkonium shown as a unitary cut across two hard gluons. The 
soft gluons (ending in a "x") are described by the field in the vacuum condensate (0|G^|0). 



The leading nonperturbative corrections to the three-gluon annihilation rate of a heavy ^Si quarko- 
nium are described by the mechanism shown in Fig|3l The soft gluon field converts the color-singlet 
^Si state into a color-octet ^Po,2 state through the El chromoelectric interaction, or into a color-octet 
^5*0 state due to the chromomagnetic Ml interaction. The final state in either of these transitions then 
decays into two hard gluons. The contribution of this mechanism can be calculated as the imaginary 
part of the graphs shown in Fig I3] given by the unitarity cut across the hard gluons, while the averaging 
of the quadratic expression in the soft gluon field makes the effect proportional to the gluon condensate 
(0|G'|0). 

It is quite essential that both the contributions of the ^Pj and ^5*0 are comparable even for a non- 
relativistic quarkonium. Indeed, the chromomagnetic Ml transition to the ^Sq state is suppressed by 
the inverse of the heavy quark mass, l/mq, in the amplitude, while the amplitude of annihilation of 
the ^Po,2 states contains an extra factor of l/mq in comparison with that for the ^5*0. Furthermore, the 
contribution of the ^Pj intermediate states is proportional to the vacuum condensate of the chromoelec- 
tric field, (0|i?^|0), while that of the ^5*0 intermediate state is given by the chromomagnetic condensate, 
(0|-B^|0). The Lorentz invariance of the vacuum however requires that these two condensates have 
opposite sign [12]: 

(015^10) = -(0|i2|0) = i(0|G2|0) . (41) 

The important feature of the discussed mechanism is that its dependence on the wave function of 
the "^Si state is not reduced to '?/'(0), but rather is determined by the corresponding overlap integrals. 
Therefore such mechanism generally breaks the similarity between the decays of different ^5*1 states. 
The same behavior however implies that the overlap integrals can be reliably evaluated only for an 
asymptotically heavy quarkonium, whose dynamics is dominated by the Coulomb-like potential due 
to the gluon exchange between the quark and the antiquark. In this limit one finds for the relative 
correction to the decay rate of an n^Si state the expression [85] 

r^^^i^ 2ii9(7r2-9)a,(mQ)A;>4 ' ^^'^ 



where a„ are numerical coefficients given by 

64 (^2 - 64) 



On = 
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(81^2- 256) (81n2-64) 
and fc" is the 'Bohr momentum' for the n-th state determined from the relation 



4 

kn = iz—rriQasikn) . (44) 

O 71/ 
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It should be mentioned that the normahzation gluon field G^i/ used here corresponds to the gluon term 
in the QCD Lagrangian of the form L = — (1/4(7^) G^, so that numerically the gluon condensate is 
(0|G2|0) ^O.TGeV^. 

As previously discussed the chromoelectric and the chromomagnetic contributions have opposite sign, 
so that the coefficients a„ in Eq. fH3|) . have the form of a difference: a„ = — a„(M). Furthermore, 

the coefficients a„(i?) at > 3 and a„(M) at n > 2 are very small. This is a consequence of a relatively 
weak Coulomb-like interaction in the octet state. If the latter interaction were neglected altogether (the 
Nc ^ oo limit) those coefficients would also vanish, and the only remaining nonzero would be ai{E) = 
a2{E) = 28 and ai(M) = 27. The real world values given by Eq.(^ are ai = 22.858 - 18.897 = 3.961 
and 02 = 9.393 — 0.323 = 9.070. These values illustrate that for the ground state the discussed 
correction is very sensitive to the details of the quarkonium dynamics, so that any extrapolation of the 
asymptotic formulas down to the realistic bottomonium and charmonium necessarily is quite uncertain. 

If such extrapolation is done by using ki ^ 1 GeV for the T resonance and ki ~ 0.5 GeV for J/ip, 
one would estimate the relative effect in the hadronic annihilation of T as about 0.5%, and about 5% 
for the J/ip resonance. If however the cancellation between the chromoelectric and chromomagnetic 
terms does not occur for these realistic quarkonia, the effect can be few times larger. 

A similar mechanism, with one of the hard gluons replaced by photon, describes the nonperturbative 
correction to the rate of the decay ^5*1 — > 7 + 2g. The relative correction in this case is three times 
smaller than in Eq. (H2l) . 

4 Radiative Transitions 

The treatment of charmonium as a nonrelativistic system suggests that one can apply the standard 
multipole expansion in electrodynamics to calculation of the transitions between charmonium levels 
with emission of a photon. The leading terms in this expansion are the El and Ml terms, which can 
be described for the discussed here transitions by the corresponding terms in the Hamiltonian: 

Hei = -Cc e{f-E) and Hmi = -f^c (A ■ 5) , (45) 

where Cc = 2/3 is the charge of the charmed quark, fi^ = Cc e/(2mc) is its magnetic moment, E and B are 
standing for the electric and the magnetic field, and A being a spin operator defined as A = cii — (T2 with 
ai and (T2 acting respectively on the quark and the antiquark. The electric dipole term is responsible 
for the transitions between the S and P states with the same spin S of the quark pair, while the Ml 
term describes the transitions between S = 1 and S = states with the same orbital momentum L. 

4.1 El Transitions 
4.1.1 tp'^-fXcj 

The rate for the transitions from a ^Si state to ^Pj induced by the El term in Eq. fHSj) is given by 

T^S, 7 'Pj) = (2 J + 1) ^ el a \Ips\' , (46) 
where uj^ is the energy of the emitted photon, and Ips is the radial overlap integral: 

POO 

Ips = {P\r\S) = / r"" Rp{r) Rs{r) dr , (47) 
Jo 

with Rs,p{r) being the normalized radial wave functions for the corresponding states. 
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Clearly, the overlap integral has dimension of length, and its particular value is somewhat model 
dependent. (With certain restrictions following from general Quantum Mechanical considerations ^S^l II].) 
One can however assess the overall validity of the nonrelativistic formula for charmonium by 
comparing with each other the experimentally known rates of the transitions from ip' to the Xcj levels 
with different J. This can be done in terms of the value of the overlap integral, \Ips\-, extracted from 
the rate of each of the transition. Proceeding in this way one finds 

{0.37fmfrom ip' Xcol 
0.39 fm from tp' ^ Xal (48) 
0.45 fm from tp' — * Xc2 7 ? 

so that the extracted values of the overlap integral agree with each other within the expected accuracy 
of the expression P6|) . and the absolute value of the radial integral also agrees well with a general 
understanding of the typical size of charmonium. 

4.1.2 Xcj^lJ/i^ 

The transitions from ^Pj levels to a ^5*1 state are described by the expression for the rate 

rCPj-.^'S,) = ^elau;'^\Isp\' . (49) 

Performing a similar to previous extraction of the overlap integral from the experimental data on the 
rates of the Xcj ^ J/'4' decays, one finds the overlap integral for the IP IS transitions as 

{0.36fm from Xco ^ J/'ip'y 
0.38 fm from Xa ^ ^7 (50) 
0.37 fm from Xc2 J/i^l ■ 

Thus one can conclude that the observed El transitions in charmonium present us with no unusual 
behavior. 



4.1.3 K^-fr]c 

The decay ^ 'yrjc is crucial for identifying the charmonium ^Pi resonance in the experiments [221 ED]- 
A straightforward application of the El transition formula results in the expression for the rate of this 
decay 

TCPi^l'So) = ^elau'^\Isp\\ (51) 

where the overlap integral is the same as is just estimated from the Xcj ^ jJ/ip decays, up to the 
relativistic corrections 0(f^/c^). Thus using the extracted value of \Isp\ one can readily predict the 
width of the decay in absolute terms: 

T{h^ ^-fVc)^ 0.65 MeV . (52) 

As discussed previously, the hadronic width of the he is quite uncertain. The logarithmic formulas 
and the data on the hadronic width of Xci indicate that hadronic decay rate of he can be as large as 
1.5 MeV, in which case the radiative decay branching ratio should be about 30%. It is believed that 
this branching ratio, with all the uncertainties involved, should be in the range between 30% and 50%. 
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4.1.4 Relativistic Effects 

In the next order of the relativistic expansion the leading amplitudes of the considered (dominantly) El 
transitions, receive 0(f^/c^) corrections from the magnetic quadrupole M2 and the electric octupole 
E3 terms [871 188] . These corrections are expected to be suppressed relative to the leading term by a 
factor f^/c^ ^ 0.2. In the total decay rate these terms do not interfere with the leading one, so that 
their contribution is only of order few percent and hardly can be measured. However these amplitudes 
generally interfere with the El term in the angular distribution. In particular for the IP IS 
transitions the M2 and £"3 contributions can be studied by measurements of the angular distributions 
and helicity amplitudes in the processes pp — > Xci,2 ^ j J/'ip- The most precise available measurements 
were performed by the E835 experiment [89j. The current averages |2 7 j for the relative values of the 
amplitudes are: M2/\E1\ = -0.13 ±0.05 and E3/\E1\ = 0.01ll|]:[i^^ from the Xc2 -^-fJ/^jj decays, and 
M2/\E1\ = -0.002l[5:[j?? from xa ^iJ/^. 

Thus the data on the radiative decay of the Xc2 generally support the expectations: the suppression 
of the M2 amplitude is indeed of order f^/c^, while the E?) amplitude was predicted to be small[8_8]. 
On the other hand the result of the extraction of M2 from the Xci decay is neither compatible with that 
for the Xc2, nor does it comply with theoretical expectations. Clearly, this situation strongly suggests 
that further studies of the relativistic terms are needed. 

4.2 Ml Transitions and the Puzzle of ijc 

The radiative transitions induced by the Ml term (Eq. fHSl) ) appear to be simpler than the El ones, 
since the Ml interaction does not contain any coordinate dependence, while the spin matrix elements 
are trivial inasmuch as the spin and coordinate degrees of freedom are factorized. Phenomenologically 
the most interesting are the Ml transitions J/ip ^ 'j rjc and ip' ~^ IVc, whose rates are given by 

4 

T(n'S, ^ 7 "^'^o) = i^ela^ \Imn\' , (53) 

3 

where Imn is the overlap integral for unit operator between the coordinate wave functions of the initial 
and the final charmonium states. In the leading nonrelativistic order, where there is no effect of 
the spin-spin interaction on the coordinate wave function, the overlap integrals are determined by 
the orthonormality condition for the coordinate wave functions: Imn = Smn- Using for an estimate 
rric = 1.4 GeV one thus calculates the rate of the transition J/ip ■jrjc as T{J/ip IVc) ~ 3.3 keV, 
while the overlap matrix element vanishes for the 25* IS transition ip' "yVc- The latter decay 
becomes possible due to the relativistic effects, but at present any estimate of the actual rate would be 
quite unreliable. 

In reality both decays have very similar rates: T{J/il) —>■ jric) ~ 1.2 keV and r(^/'' —^IVc) ~ 0.9 keV. 
Although for the latter transition the small value of the rate is not surprising and can be readily 
attributed to the suppression of the relativistic effects, the failure of the seemingly robust theoretical 
expectation for the former transition is quite paradoxical. 

One can certainly attempt to explain the discrepancy of almost a factor of 3 between the experimental 
value and the straightforward theoretical estimate for the rate of the transition from the J / in terms of a 
larger charmed quark mass, modified magnetic moment of the quark, and enhanced relativistic effects, 
or a combination of these factors. However such a 'post diction' would inevitably have to be quite 
contrived. Furthermore, Shifman has presented a phenomenological argument |90] . which eliminates the 
uncertainty related to the charmed quark mass and magnetic moment, and still leaves a considerable gap 
between the theoretical and the experimental values of the rate. The argument is based on applying 
the standard dispersion relation to the amplitude of the decay rjc — > 77. Namely, this amplitude is 
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described by just one form factor F: 

A{rj^ ^ jj) = F e^-^'- F^'J F^^J (54) 

with -Fj^'^-* being the field tensors for the two photons. The form factor is a function of the squares of 
the three 4-momenta, involved in the process: F(m^^, kf, /C2), and the width is determined by its value 
for on-shell photons: T{ric 77) = |F(m^^, 0, 0)pm^y47r. One can further use the dispersion relation 
for the form factor in one of the photon momenta: 



F{ml, 0, k') = - J j—^ lmF{ml, 0, s) , (55) 

and express the imaginary part in the integrand in terms of contribution of real intermediate states. 
These states with quantum numbers J^*" = 1 include the J/ip resonance, higher charmonium reso- 
nances: '?/'(3770) and so on, and the continuum of the states with charmed D{D*) mesons. 

The contribution of the J/ip is expressed in terms of the product of the amplitudes A(?7c 
7 J/ip) A{J/ip 7*(s)), where the coupling of J/ip to a virtual photon is obviously related to its decay 
width into e~^e~ . The contribution of the higher charmonium resonances is vanishing, as discussed, in 
the leading nonrelativistic approximation and is indeed very small phenomenologically. Namely, the 
effective overlap integral J21, that would fit the observed rate of the decay ip' IVc-, is numerically 
only about 0.04. Similarly the contribution of the continuum to the integral in Eq. (l55|) can also be 
argued[_90j to be quite small compared to the contribution of the J/ip. Therefore the dispersion relation 
( !55ll should be saturated to a good accuracy by the J /ifj pole thus relating the amplitudes of the decays 
Vc — 77, J/i^ ^ e~^e~ and J/ip —>■ 'yrjc and the masses of rjc and J/V', i-e. without the need to invoke 
the mass of the charmed quark. The resulting relation between the rates reads as 

r(T/i \ 2 T{r],^ 77) ^J/^b^-. "^l V n ^ M (r^a^ 

r(j/^ - 7^/.) = - ^^j^^ ^ + ^''^ 

Using the experimental data in the r.h.s. of this relation, one finds the central value of the discussed 
rate (modulo the QCD corrections) as T{J/ip •yrjc) = 2.9 keV, which is pretty close to the previous 
simple nonrelativistic estimate, and is still very far away from the experimental value of this rate. 

It may well be that the puzzles of the 7]^. its large total width, small branching fraction of decay into 
77 and the suppression of the decay J/ip ^ jrjc are all tied together and possibly can be resolved by 
a strong mixing in the ^ channel [5 9]. so that the rjc has a sizable admixture of light quark and gluon 
states. Qualitatively such admixture would enhance the hadronic decay rate and suppress the radiative 
transition from a pure charmonium state J/ip. In terms of the Shifman's argument, there would arise 
an additional contribution to the imaginary part of the form factor F from 1 states lighter than J/ip, 
which would generally violate the relation fl56l) . However a quantitative description of these phenomena 
cannot be offered at present. 



5 Hadronic Transitions Between Charmonium Resonances 

Similarly to the radiative transitions between heavy quarkonium levels, caused by the interaction of 
quarks with photons, the hadronic transitions arise through the interaction of the heavy quarks with 
gluons and the gluons materializing as light mesons, the pions and rj. Also in complete analogy with the 
radiative transitions, the interaction of a nonrelativistic quarkonium with the gluon field in hadronic 
transitions can be described within the multipole expansion in QCD[91J. The leading terms, that 
will be important in our further discussion of the realistic processes are the chromoelectric and the 
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chromomagnetic dipoles, El and Ml and the chromomagnetic quadrupole M2. The corresponding 
terms in the Hamihonian can be written as 

Hei = —Cr-E'^iO) , Hmi = -^Ci^-B^ , and Hm2 = -{^mg)-' C {D,B,{0)y , (57) 

where C"' = — is the difference of the color generators acting on the quark and antiquark (e.g. 
ti = A"/2 with A" being the Gell-Mann matrices), r is the vector for relative position of the quark and 
the antiquark, 5* = {aq + aQ) /2 is the operator of the total spin of the quark- antiquark pair, and D is the 
QCD covariant derivative. Finally E and B are the chromoelectric and chromomagnetic components 
of the gluon field strength tensor. 

An important difference of the hadronic transitions from the radiative ones is that the physical 
amplitudes arise in at least the second order in the interaction with the gluon field due to the color 
requirements. Considering quarkonium as a compact object interacting with soft gluonic fields, one 
can further approximate the quarkonium transition in the second order in the interaction Hamiltonian 
in terms of a local colorless glounic operator, which operator produces the light mesons in the actual 
hadronic transition. Thus the full transition amplitude factorizes into the heavy quarkonium part 
determined by the terms in the multipole expansion (Eq. fl57p ) and the production amplitudes of light 
mesons by the gluonic operators. The heavy quarkonium part is to a great extent model dependent, and 
here we discuss this part only by relying on the general dynamical properties, while the production of 
soft light mesons by local gluonic operators can be described using chiral algebra and certain low-energy 
theorems in QCD. In what follows we first concentrate on the latter description, and then apply the 
results to specific hadronic transitions. 



5.1 Production of Light Mesons by Gluonic Operators 

The structures arising in the second order in the interaction terms in Eq. (!57l) are quadratic in the gluon 
field strength tensor. Therefore the interesting for the present consideration amplitudes are those for the 
production of one or two light pseudoscalar mesons by the operator of the general form C^^C^^. As will 
be described, in the low-energy limit the two-meson production by this operator is determined [92| [93] . up 
to a small constant, from the chiral algebra and the QCD anomaly in the trace of the energy-momentum 
tensor [9^ [95l [96l [97] . while the one-meson production amplitude is fully fixed by the anomaly in the 
light quark axial current [HSl ES]- Furthermore, the matrix element for one- meson production by the 
operator containing one extra covariant derivative, G^yDG\rj is entirely determined |92[ llUUj (in the 
same low-energy limit) by the fixed matrix element without the derivative, and is therefore relatively 
well understood. The knowledge of the production amplitudes for the light mesons in fact results in 
relations between the observed transitions which appear to hold quite well when compared with the 
data. 



5.1.1 Two-pion Production Amplitude 

Let us consider first the two-pion production amplitude {tc^ iPi)Tc^ iP2)\G"^^G1^\0)^. In the leading chiral 



limit the momenta pi and p2 of the pions as well as the pion mass are to be considered as small 
parameters, and the expression for the amplitude, quadratic in these parameters, can be written in the 
following general form 



- {^^iPi)^-iP2)\G%GlJ0) = \X {p, .p,) + Y {pI + pI 



i9^l\9ua - gfiagux) + Zt^yXa , (58) 



•^The charged pions are considered for dcfiniteness. The amphtude for the neutral pion pair production is trivially 
related by the isospin symmetry. 
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where the structure 



- {plfiP2a + PlaP2fi) QuX " {PluP2\ + Pl\P2u) Qfia " {Pl ' P2) (fi'/.Afi'i/a " gfiaQux) (59) 

has zero overall trace: tf^^'^" = 0, and X, Y, and Z are yet undetermined coefficients. The form of 
the amplitude in Eq.( l58i) is uniquely determined by the symmetry (with respect to the indices) of the 
operator G^^Gx^j and by the Adler zero condition, which requires that the amplitude goes to zero when 
either one of the two pion momenta is set to zero and the other one is on the mass shell. One can also 
notice that the proper index symmetry and the Adler zero condition also automatically ensure that the 
amplitude is C even, i.e. symmetric under the permutation of the pion momenta: pi ^ p2- 

The coefficients X and Y in Eq. flSHl) are in fact determined [92] by the anomaly in the trace of the 
energy-momentum tensor in QCD. Indeed, in the low-energy limit in QCD, i.e. in QCD with three 
light quarks, one finds 

0', = G^^G^^-^ + Y: m,{qq) , (60) 

q=u,d,s 

where 6 = 9 is the ffist coefficient in the beta function for QCD with three quark flavors. The ffist term in 
Eq. fl60l) represents the conformal anomaly, while the quark mass term arises due to the explicit breaking 
of the scale invariance by the quark masses. On the other hand, the matrix element of the energy- 
momentum tensor 6^^ over the pions: 9^^{pi,p2) = {t^^ (pi)'^^ {P2)\(^fiu\0) is fully determined [5^1^ [TUT] 
in the quadratic in pi,p2 and m^r order by the conditions of symmetry in /i and z/, conservation on the 
mass shell {{pi + P2Y 6^y{pi,p2) = Q aX p\= p'\ = rri^), normalization {O^uip, —p) = '^P^lPu at p^ = m^), 
and the Adler zero condition {O^uij)^ ^)\p'^=ml = 0): 

0^u{Pl,P2) = [{Pl ■ P2) +pI+pI- ml\ g^,^ - Pi^P2u - PlyP2t, ■ (61) 

The equations (135]) and flMl) allow for the pion momenta to be off-shell in order to demonstrate the 
Adler zero. However in what follows only the amplitudes with pions on the mass shell will be considered, 
so that it will be henceforth implied that pi = pi = to^. In particular one finds for the trace of the 
expression in Eq. fl61l) 

e^^{p,,P2) = 2{p,-p2)+Aml. (62) 

Furthermore, the quark mass term in Eq. (l60l) . corresponding to the explicit breaking of the chiral 
symmetry in QCD corresponds to the same symmetry breaking by the pion mass term in the pion 
theory. Thus one finds to the quadratic order in m^: 

(tt+tt^I ^ mq{qq)\Q) = ml , (63) 

q=u,d 

while the term with the strange quark makes no contribution to the discussed amplitude. 

Combining the formula in Eq. fl60l) for 6'^ with the expressions (162]) and (163]) one finds the matrix 
element of the gluonic operator over the pions in the forno 

- {n^{pi)n-{p2)\G';,,G'^nO) = [2 (pi ■ P2) + ^ ml] (64) 

which thus allows to determine the coefficients X and Y in Eq. (IS5]) : X = 167r^/(36) and Y = 3X/2 = 
87rV6. 



"'it can be mentioned that this relation, taking into account the pion mass, was used in Ref. [102] . and was also derived 
in a particular chiral model in Refs. [1031 1104] . 
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The coefficient Z of the traceless part in Eq. fl58|) cannot be found from the trace relation fl60|l . 
Novikov and Shifman |93j estimated this coefficient by relating this part to the matrix element of the 



traceless (twist-two) energy-momentum tensor of the gluons only: 6*^^ = Aixds [—G°a\G°^^\9iJiv G^^G 



a Xo 



Zt.xu'' = 47ra, (7r+(pi)7r-(p2)|e|0) 



(65) 



They then assume that the matrix element of the twist-two operator is proportional to the traceless 
part of the phenomenological energy momentum tensor of the pions, 



"(P2)ieio) 



PG 



{Pl ■ P2) Qpu - PlpP2u - PluP2p 



(66) 



with the proportionality coefficient interpreted, similarly to the deep inelastic scattering, as "the fraction 
of the pion momentum carried by gluons". They further introduce a related parameter k = basPc/iQ'^) 
and conjecture that numerically k ^ 0.15 — 0.20. However, for the purpose of considering the two-pion 
transitions it is not necessary to pursue the interpretation of k as being related to the gluon structure 
function of pion, but rather it can be treated as a phenomenological parameter which can be determined 
from the data. 

Summarizing the results so far in this section one can write the expression for the general matrix 
element (1581) for on-shell pions as 



(vr+(pi)7r-(p2)|F;,Ft|0) 



3b 



(g^ + m^) (g^xg. 



Qpugux) 



9 

2 ^ tpuXa 



(67) 



where q = pi + p2 is the total four-momentum of the dipion. 

Few remarks are due regarding effects of higher order in as- The trace term in Eg. (1671) receives no 
renormalization, provided that the coefficient b is replaced by ^{as)/^^ with P{as) = ba^ + O^al) being 
the full beta function in QCD. This modification however only affects the overall normalization of the 
trace part, and can in fact be absorbed into the definition of the heavy quarkonium amplitudes. On the 
contrary, the relative coefficient of the traceless term in Eq. fl67|) . i.e. the parameter k, does depend on 
the normalization scale, which scale is appropriate to be chosen as the characteristic size of the heavy 
quarkonium |93]. However, given other uncertainties in the analysis of the two-pion transitions, the 
slow logarithmic variation of k is a small effect. 



5.1.2 Production of rj by Gluonic Operators 

The amplitude of production of the rj meson by the quadratic gluonic operator is described by only 
one form factor: {rj\G ^yG x„\0) = Atf^^xa and is therefore entirely determined by the expression [OSl ES] 
following from the chiral algebra and the anomaly in the divergence of the singlet axial current in QCD, 

e""'^ (r/IG^.GA.lO) = 167r2 ^F,mJ , (68) 

where is the r] 'decay constant', equal to the pion decay constant ~ 130 MeV in the limit of exact 
flavor SU(3) symmetry, and is likely to be larger due to effects of the SU(3) violation. 

A more interesting case is presented by the matrix elements relevant to discussion of transitions in 
quarkonium involving the M2 interaction from Eq. (l57|l and containing an extra covariant derivative: 
{ri{p)\Gfj_yDpGx„\0) and {rjl^DpG^^YGI^lO). It turns out that these matrix elements are also entirely 
determined by the relation (I68p . The possibility of the reduction of the structures with extra derivative 
to the expression in Eq. (1681) follows from the general theory [TU5] . The reasoning in this specific case [100j 
makes use of the following identity, valid for arbitrary four-vector p: 

Pp^p,uXa PX^fjLPpa PcT^p,upX P^^upXa ~\~ Pu^ppXa i (69) 
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where the convention €0123 = 1 is assumed. The antisymmetry of the field tensor G^i/ then allows 
one to write the general form of the first of the discussed matrix elements in the linear order in the t] 
momentum p in terms of two scalars X and Y: 

) • (70) 

The third structure, allowed by the symmetry and proportional to {p^e^pXa — Pu^fipXa), is reduced to 
the first two due to the identity flB^ . Furthermore, applying in Eq. flTD]) the equations of motion (the 
Jacobi identity): DpGxa + D„Gp\ + DxG^p = 0, one arrives at the relation X = 2Y. 

Likewise, writing the second of the discussed matrix elements in terms of two scalars X and Y as 

z{ri{p)\{DpG,,rGUO)=X 

apu Pu^Xapp ) , (71) 

and applying the Jacobi identity, one finds the relation X = 2Y. 

The sum of the expressions (!70l) and (17T]) should combine into a total derivative, i.e. the sum should 
be proportional to Pp. This is possible due to the identity (!69|) under the condition that Y = Y, so that 
all the considered scalar form factors are expressed in terms of one of them, e.g. in terms of X-^ 

X = X, Y = Y=^X. (72) 

Using this relation and contracting the sum of the expressions f lTOj) and fl7T]) with ^e^^^^^ the form factor 
X is identified from the equation (|68|) as 

^ = e"""- {v\Gp.Gx.\0) = ^^v< ■ (73) 

5.1.3 Production of Single tt^ by Gluonic Operators 

The production of a single neutral pion by gluonic operators obviously requires a breaking of the isotopic 
symmetry. Within the chiral approach the isospin breaking originates in the different masses of the up 
and down quarks, mu 7^ w^d. Allowing for this mass difference, the gluonic matrix elements for a single 
vr" can be considered in essentially the same way as those for the 77 meson. In particular the 'master 
equation' for the anomalous amplitude in the case of pion takes the form [98] 

e'^^^'^ {ri\Gp.GxM = IGvr^ " ^" . (74) 

rrid + vfiu 

Thus the simple rule of conversion from the amplitudes of the processes with 77 emission to similar 
processes with single 7r° can be done by replacing the coefficient — > (m^— m„) / {nid+niu) 

and replacing the r] momentum with that of the pion, ^ p-k- Naturally, these replacement rules can 
be also viewed in terms of a fixed rj — mixing. 

5.2 Two-pion Transitions 
5.2.1 The Transition ip' 7111 J ftp 

The general soft-pion relations of the chiral algebra require that the amplitude of a two-pion transition 
is bihnear in the components of the four-momenta of the pions p!06| I107j in the chiral limit. For the 

^An alternative derivation of two of these relations, namely X ^ X and Y = Y, would be by arguing that in the 
particular amplitudes ([70)) and ([7T|) the operators G° and [DG)"" can be considered as commuting with each other, so 
that the expressions ([70|) and ([7T|) differ only by re-labeling the indices. 
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most studied transition %p' — ^ TrnJ/tp this implies that in the hmit of soft pions the amphtude can be 
generally parametrized as |lU6j 

A{^' -^TT+TT-J/iP) = [A{q^-2ml) + Xml\{e,-e2) + BE,E2{e,-e2) 

+ C [(pi-ei)(p2-e2) + (p2-ei)(pi-e2)] , (75) 

where e^* and €2 are the polarization amplitudes of the initial and the final vector resonances, pi and 
P2 are the 4-momenta of the two pions, Ei and E2 are their energies in the rest frame of the initial 
state, and q = pi + p2 is the total 4-momentum of the dipion, so that = m^^. Finally, A, B, C, and 
A are the form factors, which should be considered constant in the soft pion limit, but are generally 
functions of kinematic variables beyond this limit. The spin-dependent term, proportional to C should 
be suppressed inasmuch as the charmed quark can be considered as heavy, since the spin-dependent 
interaction is proportional to the inverse of the heavy quark mass. The remaining constants in Eq. fl75l) 
can be related to the parameters in the leading order of the multipole expansion in QCD. 

The two-pion transition between ^Si states is generated in the second order in the leading El term 
(Eq. fl571) of the multipole expansion, so that the amplitude of the process can be written as 

A{^' ^ yr+TT- J/^) = ^-{n+n-\E^E^\0) ajf , (76) 

(12) 

where a-^ can be termed, in complete analogy with the atomic properties in electric field, as the 
transitional chromo-polarizability of the quarkonium. In other words, the ip2 ^1 transition in the 
chromo-electric field is described by the effective Hamiltonian 

Heff = -^alf E^E^ , (77) 
with the chromo-polarizability given by 



cylf = ^{lS\engr,e\2S) , (75 



where Q is the Green's function of the heavy quark pair in a color octet state. The latter function is 
not well understood presently, so that an ab initio calculation of the chromo-polarizability would be at 
least highly model dependent. Generally aij is a symmetric tensor. In the leading nonrelativistic order 
the coordinate and spin degrees of freedom factorize, so that the discussed chromo-polarizability for a 
transition between S wave states is actually reduced to a scalar: 



(12) 

alj = a 



^''^ 5,, (6-1 ■ e2) . (79) 



The matrix element for production of the pion pair by the quadratic chromoelectric operator is then 
found from the general expression fIBTl) . and the result for the amplitude of the decay ip' tt^ti^J/'iIj 
ispiigglfTUg] 



A{iP' n+n^J/iP) = (80) 

pfi pu' 



(g + m. 




ei ■ €2 



where P is the 4-momentum of the initial tp' resonance, and the tensor corresponds to a D wave 
motion in the cm. frame of the dipion and is written in terms of the relative momentum r = pi — p2 
of the two pions as 

1 / 4m^ \ 

^f,u = r^r^ + - ( 1 - — ^ j (g2 g^^ - q^q^) . (81) 
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Figure 4: The spectrum of the dipion invariant masses m7r,r in the decay ip' tt+tt J/ip as 
given by the amphtude (IHOl) with k = 0.2. 



The constants in the parametrization fl75l) can be thus found as 



A 



An' 



a 



(12) 



^ 3 

1 + -K 

4 



B = 6k 



■ a 



(12) 



A 



12k 
4 + 3k 



f82) 



while the parameter C is clearly equal to zero, since the El interaction carries no spin dependence. 
It can be also noted that Eq. flSU]) differs from (1751) at C = in that it contains only two parameters, 
while the third, A is fully fixed as shown in Eq. fl82l) . This is due to the chiral constraints and the 
conformal anomaly (Eq. fl60l) ) fully determining the terms of order relative to those quadratic in the 
pion momenta. 

One can readily see from Eq. fl80l) that the overall rate of the decay is determined by the chromo- 
polarizability a^^'^\ while the shape of the spectrum depends on the parameter k. Furthermore, the last 
term in the braces in Eq. flHOj) is proportional to k and describes the D wave motion in the cm. frame 
of the two pions, which D wave correlates with the motion of the dipion as a whole in the rest frame 
of the decaying ijj'. Thus the formula flHOj) based on the leading El term of the multipole expansion 
predicts [nS] a relation between the amplitude of the D wave and the parameter in the shape of the 
spectrum in the decay. The most detailed to date experimental study of the decay ip' n^n^J/ip 
using this formula was done by BES[109]. The fit to the spectrum of the invariant mass of the produced 
dipion resulted in the value k = 0.186 ± 0.003 ± 0.006, while the fit to the ratio of the D and S wave 
amplitudes from the angular distribution gave k = 0.210 ± 0.027 ± 0.042. Clearly, the consistency of 
these two values implies that the discussed approach correctly predicts the ratio of the D wave in terms 
of the sub-dominant term proportional to k in the S wave amplitude. It can be also noted that at 
K ^ 0.2 the contribution of the D wave to the total rate is quite small - only about 2% |109j . The 
familiar characteristic shape of the spectrum of the two-pion invariant masses is shown at k = 0.2 
in FigjH Comparing the integral over this spectrum with the experimental total decay rate one can 
evaluate the transitional chromo-polarizability as a'-^^-' ~ 2 GeV~'^. 

It can be also noted that in the discussed picture the pions are produced in pure isoscalar 7 = 
state, which implies fixed relation between the differential rates of emission of pairs of neutral and of 
charged pions dTlip' — >■ ir^n'^J/ip) = 0.5dr{ilj' — ^ ir'^n' J/ip). The most precise test of this relation has 
been done by the CLEO-c experiment [llOj . which resulted in the measured ratio of the total decay 
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rates 0.4924 ± 0.0047 ± 0.0086. 



5.2.2 Effects of the Final State Interaction Between Pions 

So far the amplitude of the two-pion transition was considered in the chiral hmit. The formula in 
Eq.( l67l) is exact in the leading chiral order, i.e. as far as only the quadratic terms in the pion momenta 
and mass are concerned. In particular this expression receives no corrections due to the final state 
interaction (FSI) between the pions. The latter interaction however can give rise to the terms whose 
expansion starts with the fourth power of momenta and the pion mass, and generally can modify the 
amplitude at momenta of the pions relevant for actual transitions in quarkonium. The effects of FSI 
in chiral treatment of the two-pion transitions in quarkonium were a matter of concern ever since the 
earlier theoretical analyses |106j . The effect in the phases of the amplitudes is well known: these phases 
for the production amplitudes are equal to the two-pion scattering phases in the corresponding partial 
waves: S = \S\ exp(z5o), D =\D\ exp{i62), where the / = phases for the 5* wave, 5o, and for the D 
wave, §2 are quite well studiedjj It is also generally estimated both on theoretical and phenomenological 
grounds that the FSI corrections are not big (at most 20 - 25%) in the transitions ip' — > uttJ/iIj and 
T' ttttT. (For a discussion see the review Some phenomenological arguments in favor of such 

estimate will also be discussed here after presenting a theoretical estimate of the onset of the higher 
term in the chiral expansion. 

The interaction of pions at low energy in the D wave is quite weak, so that any modification by FSI 
of the D wave production amplitude can be safely neglected, and only the modification of the S wave 
amplitude is of interest for present phenomenology. The imaginary part of the correction at > 4m^ 
is found from the unitarity relation in terms of the isospin 1=0 tttt tttt scattering amplitude T[q'^) in 
the wave as 

,„(,S)^XrMl|!)5. (83) 
V q'^ lovr 

The amplitude T{q^) is well known in the chiral limit, i.e. in the quadratic in q and approximation, 
since the work of Weinberg [112] . In the normalization used here this amplitude has the form 

T(g^) = ^-^j^ , (84) 

where ~ 130 MeV is the 7r+ fi^u decay constant. Clearly, the expression in Eq. (l83l) is of the 
fourth power in q and m. 

The real part of the correction to the 5* wave production amplitude 6S can then be estimated from 
Eq. fl83l) using the dispersion relation in for the amplitude S. In doing so one should set the condition 
for the subtraction constants that this real part does not contain quadratic (and certainly also constant) 
terms in q and m, since these are given by Eq. fl^7|) . After these subtractions the remaining dispersion 
integral is still logarithmically divergent and contains the well known 'chiral logarithm', depending on 
the ultraviolet cutoff A, which is usually set at A ~ 1 GeV, i.e. the scale where any chiral expansion 
certainly breaks down. Using this approach one can estimate |108j the corrections in the S wave part of 
the amplitude given by Eq. (!80|) as 

A{tlj' 71+71- J/tP) = (85) 



2 



{q' + ml){l + ^,)-K 1 




^It can be mentioned that the analysis [109j of the data on the -0' tt+tt^ J/i/j decay does not take into account the 
relative phase between the S and D wave pion production amplitudes. Thus it would be interesting to know whether 
including the phase factor in the angular analysis, produces a significant impact on the results. 
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where the correction terms .^i and ,^2 are given by 



167r2 F2 (g2 + ^2 ) ' IQtx V ' l^'^^ 

and 

167r2F2(g2 + 2m2) ml IGti V ' 

where the non-logarithmic imaginary part is retained for reference regarding the normahzation. The 
lower limit under the logarithm is generally a function of both and m^, however any difference of 
this function from the value used in Eqs.f l86p and flHTl) is a non- logarithmic term, i.e. beyond the 
accuracy of these equations. Since ml is the smallest of the two parameters in the physical region of 
pion production, it can be expected that using this parameter provides a conservative estimate of the 
effect of FSI. 

Estimating the corrections in Eq. fl86|) and Eq. fl87|) . one finds that at the lower end of the physical 
phase space, i.e. near = 4m^, these terms do not exceed few percent. Thus the corrections only 
weakly modify the normalization of the pion production amplitude near the threshold. A theoretical 
extrapolation to higher values of is problematic, and, most likely, one would have to resort to using 
actual data on the dipion spectra in order to judge upon the significance of deviation from the essentially 
linear in behavior of the amplitude described by Eq. (l80|l . A quantitative estimate of the deviation 
from this behavior has been attempted [lllj using the data |113j on T' — >■ -7r"'"7r~T by parametrizing 
the deviation as a factor (1 + g^/M^) in the amplitude with M being a parameter. The thus obtained 
lower limit on M is 1 GeV at 90% C.L. 

Another phenomenological argument in favor of a relatively moderate FSI effect in the absolute value 
of the dominant S wave in the decay ip' — > t[^ti~ J /if) stems from the previously mentioned agreement 
of the observed [lOQj value of the ratio D/S with the parameter k entering the expression for the 5* 
wave and extracted from the two-pion invariant mass spectrum. Clearly such an agreement would be 
ruined if there was a significant enhancement of the S wave by FSI. 

The treatment of the FSI effects in the two-pion transition amplitude in fact reveals a certain 
inadequacy of the parametrization in Eq. (l75l) . Indeed the factors B and C are each contributed by both 
the S and D wave motion in the cm. frame of the two pions, and those contributions are differently 
modified by the FSI effects. In this situation it is more reasonable to use the parametrization jl08[ 1114] 
in terms of separate partial waves: 

p^l pu 

A{^' ^ TT+TT- J/^) = ^ (ei ■ es) + D, i^, (ei ■ €2) + q, g. e^^ + ^^.u e^" , (88) 

where e^'^ stands for the spin-2 tensor made from the polarization amplitudes of the ^5*1 resonances 



e^'^ = e^^ + e^e^ + - (ei • e^) -— - g,^ . (89) 



/ p^ 

The terms in the expression fl88|) describe an S wave and three possible types of D-wave motion: the 
term with Di corresponds to a D wave in the cm. system of the two pions correlated with the overall 
motion of the dipion in the rest frame of the initial state, the D2 term describes the D-wave motion 
of the dipion as a whole, correlated with the spins of the quarkonium resonances, and finally, the 
term corresponds to the correlation between the spins and the D-wave motion in the cm. frame of the 
dipion. Clearly, the factors S and D2 are modified by the FSI in the S'-wave dipion, while the factors 
Di and receive a modification from the two-pion D wave FSI. The partial-wave form factors are 
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expressed in terms of the parameters in Eq. flTSl) as 



S 



B 



2ml) + Xml + — B--C 
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3 go - (^0 - g^ 



1 - 



(90) 



Naturally, in the presented treatment, based on the multipole expansion in the leading nonrelativistic 
order, only the spin-independent amplitudes 5* and Di are present, as described by Eq. llHOl) and Eq. (!85ll . 



5.3 Single rj and Single Pion Transitions 
5.3.1 ij'^riJ/^ 

The decay ip' V^/ 4' arises due to the interference of the El and M2 terms in the multipole expansion 
(Eq.( l57I) ). The amplitude of the transition can thus be written as 



ijk 1 



where 



1 



(91) 



(92) 



is the quarkonium transition amplitude. In the approximation of factorized coordinate and spin degrees 
of freedom the matrix element of the spin operator is expressed in terms of the polarization amplitudes 
e*! and e2 of the initial and final quarkonium resonances, while the remaining coordinate overlap is the 
same as in the expression (178!) for the chromo-polarizability, so that the amplitude is expressed in terms 
of a*^-^^) as 



a 



(12) 



4mQ 



(93) 



The amplitude of the ri production by the gluonic operator in this expression is readily found from the 
relations (170|) - (!73l) . and is given by 



i (r/|Ef (A5fe)'^ + (AA)«Ef|0) 
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with p being the momentum of the rj meson. As a result one finds the following simple formula for the 
amplitude of the 7] transition 
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Comparing the latter formula with Eq. flHOl) one finds that the quarkonium part of the amplitude, 
a^^'^\ cancels in the ratio of the amplitudes for the t] and the two-pion emission and the ratio is essentially 
determined by the two anomalies in QCD: the conformal anomaly and the one in the divergence of the 
singlet axial current [92] . This remarkable relation can be written in terms of the decay rates as 
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where the factor T describes the deviation of the two-pion transition amplitude from the limit of 
dominance of the anomaly contribution (/t 0) and of massless pion: 
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where the last term describes the small contribution of the D wave. Using k, = 0.2 and integrating over 
the phase space of the two-pion transition one can arrive at the estimate of the ratio of the total rates: 

n^'-yJ/^) = 09 ( V ^IdGeVX^ 

r(^' ^ TT+TT- J/^) USOMeVy \ m, J ' ^ ' 

Given all the uncertainties, this estimate agrees very well with experimental value 0.097 ±0.003 for this 
ratio. 

It can be mentioned as a sidenote, that being applied to the transitions in bottomonium between 
T' and T, Eq. (pB]l . using the appropriate for these transitions value k ~ 0.15, predicts the ratio of the 
rates 

r(T' ^ vr+vr-T) " i^'^ ) [uOMeV ) l^^^j ' 

corresponding to B{T' — ^ r/T) ^ 4.3 x 10~^. This prediction can be compared with the recent prelimi- 
nary data from CLEOfTTB]: B{T' ^ r/ T) = (2.5 ± 0.7 ± 0.5) x lO"! 

Another sidenote regarding the decay tp' rj J/ ip is that recently this process has been used not for 
its own sake, but rather as a precision source of the 7] mesons for high accuracy measurements of the t] 
mass [116] and its decays pTTj . 

5.3.2 ip'^TT^J/ip 

The decay ip' tt^ J/ip requires breaking of the isotopic symmetry. If the relevant to this process 
breaking is due to the mass difference between the u and d quarks the amplitude of this decay can be 
found from the amplitude of the i] transition by applying the previously discussed conversion factor, so 
that the ratio of the rates for these two transitions can be estimated ^JLSj as 

T{iIj' ^ TT^ J/ijj) _ ^ fmd-mu\^ ^^^^^ 



T{iIj' ^ 7] J/i/j) \md + muJ F^ 



In the limit of the flavor SU(3) symmetry one has F^^ = F^j. In reality it is known from comparison of 
Ft, and Fx that the presence of heavier strange quarks increases the constant F, so that is expected 
to be larger than F^. Therefore the limit F^^ = can be used as an upper bound on the ratio of the 
rates in Eq. fllOOp . The ratio of the masses of the u and d quarks, describing the explicit breaking of 
the chiral symmetry and the isospin violation in this breaking was studied years ago in great detail 
by Gasser and Leutwyler |119j . The largest value for the ratio (m^ — my)/{md + m„) allowed by that 
study is approximately 0.35. Thus the theoretical upper bound for the ratio of the decay rates in 
Eq.f llOOp is approximately 2.3%, which is still by more than 4(T below the experimental result [Hn]: 
(4.1 ±0.4 ±0.1)%. It can be also mentioned in connection with the light quark mass ratio that the well 
known Weinberg's formula [120j gives 

= = 0.285 , (101) 

vfid + ruu mio 



and results in a still smaller ratio of the decay rates if used in Eq. fllOOl) . It is certainly understood jllS 



that the formula fllOOp may receive unexpectedly large corrections from the effects of the SU(3) viola- 
tion, however such corrections can also significantly affect the analysis of the chiral phenomenology in 
Ref. |119j . and the whole subject then would have to be revisited anew. It should be mentioned that 
the largest theoretical estimate of the discussed ratio of the transition rates found in the literature |121] 
corresponds to 3.4%, which is also significantly lower than the experimental number. However, the 
latter estimate does not fully take into account the proper QCD structure of the relevant amplitude for 
the meson production by the gluonic operator. 
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It therefore looks like the isospin violation by the light quark masses is not sufficient to describe the 
data [110] and one has to assume that at least one of the charmonium states in the transition is in fact 
not a pure isoscalar, but contains a small admixture of an isovector 1 = 1 four-quark state. We shall 
further discuss such possibility in the Section 6.1.2. 



5.3.3 ip' ^Ti^hc andhc^TT^ J/^ 

The transitions between the ^5*1 and ^Pi states with emission of tt^ arise in the discussed picture through 
the interference of the El and Ml interaction terms in Eq. flS?!) . The pion is emitted in the S wave, 
and the amplitude for the transition can be written as 

^(3^1 ^vrO^p,) = {^^\E'^Bl\Q)It^kt2^ , (102) 

where I is the quarkonium radial overlap integral 

^ {P\Cirgs + gpr)C\S) (103) 
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with Qs and Qs being the octet-state Green's function in the corresponding partial waves. 

The gluonic matrix element for the single pion production is readily found from Eq. fTMl) . and one 
finds the following expression for the decay rate [122] 

res, ^ = |/pPi . (104) 

\ 3 md + rriu J n 

A numerical evaluation of the decay rate for the transition ip' —>■ greatly suffers from a poor 
knowledge of the overlap integral in Eq. fll03p . If, for an order-of- magnitude estimate, one uses typical 
r as in the radiative transitions, r ~ 0.4 fm, and also approximates Q ~ lGeV~^, then one very 
approximately estimates |/| ~ 0.15 GeV"^, and also estimates the rate from Eq. (11041) as T{iIj' —>■ 7r° he) ~ 
15 eV, which corresponds to the branching fraction of only about 5 x 10"^. The latter number is by an 
order of magnitude smaller than the measured |30j combined branching fraction B{ip' 7r° he) x B{hc 
IVc) = (4.0 ± 0.8 ± 0.7) X 10~^. It is not clear whether the estimate of the overlap / is totally off the 
mark, or the apparently enhanced decay ip' — > 7r° he also proceeds due to a small four-quark admixture 
in the ip'. 

A similar estimate of at most few tens eV is applicable to the transition from the he resonance, 
he —>■ vr° J/tp. The experimental status of this decay is still not clear. An early observation of this decay 
in the E760 experiment [28] has not been conffimed by the E835 datap^ a decade later. Clearly, an 
additional experimental input on this decay would be of a great interest. 



5.4 Chromo-polarizability and Slow Charmonium in Matter 
5.4.1 Diagonal Chromo-polarizability of J/ip and ip' 

As discussed, the hadronic transitions from ip' to J/ip are determined by the transitional chromo- 
polarizability a^^'^\ defined in Eq.flTHl). Similar diagonal quantities a*-^^-* = and a*-^^-* = a^/ determine 
the interaction of the J/ip and ip' resonances with soft gluonic fields and are of a great importance for 
description of a number of processes including the interaction of the charmonium resonances with 
nuclear matter. 

At present these characteristics of the charmonium resonances are not well known. An early 
calculation[123] of the chromo-polarizability of quarkonium ground state relied on the treatment of 
the system as being purely Coulomb-like, which is a valid limit for very heavy quarkonium, but which 
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hardly can be applied to charmonium. One guidance for the value of the discussed parameters is set by 
the transitional chromo-polarizability a*^^^^ ~ 2 GeV~'^. Namely the diagonal terms and a^i should 
satisfy the Schwartz inequality 

and they both should be real and positive since there are no intermediate states that would enter the 
second-order correlator in the El interaction over either of the vector resonances. Thus it is reasonable 
to consider the value of a''^^-' as a 'reference' benchmark for either of the diagonal terms. 



5.4.2 The Decay J/^ ^ 7r7r£+£- 

The chromo-polarizability of the J /if) resonance can in fact be measured experimentally in the decay 
J/ip ^ TT+TT" £^C.~ with a soft pion pair [124]. Using the same description based on the QCD multipole 
expansion as for the two-pion transitions, one can write the amplitude of the decay in the form 

2 m[nb) - mjj^ ^ go 

where the sum goes over the discrete states as well as the continuum, and qo = {q- P)/'mji^. In writing 
this expression it is taken into account that the soft pion approximation is only valid at go ^ ^j/i/), 
so that any recoil of the heavy quarkonium upon emission of the pion pair can be and is neglected. In 
this limit the relation between go and the total momentum / of the lepton pair can also be written as 



m 



In the chiral limit the first term (with n = 1) in the sum in Eq. fll06p dominates for soft pions, due to 
its singular behavior as 1/go. It can be noticed however that the decay amplitude itself is not singular 
due to the (even faster) vanishing of the pion production amplitude fl67|) in the limit of soft pions. In 
the 'real life' the minimal practical energy go is not much less than the spacing of the quarkonium levels, 
and the contribution of higher states mixes into the amplitude. In what follows we first consider the 
contribution of only the first term of the sum in Eq. (11061) and then discuss the effect of the higher terms. 
Keeping only the contribution of the first term in the sum in Eq.f ll06p . one can write the differential 
rate of the discussed decay in terms of the chromo-polarizability and the leptonic width TedJ/'ip) in 
the form 




dT{J/^ ^ 7r+7r-£+r) = \a^\' ^ ql - q'T,,{J/i;) dq' dqo , (107) 



where the factor is the same as defined in Eq.([9j 

In order to assess the feasibility of observing the discussed decays it can be noted that at a given 
constraint on the maximal value of go: go < A (or equivalently at a lower cutoff on the invariant mass 
of the lepton pair) the probability described by Eq. (11071) strongly peaks near the highest values of both 
g^ and go, i.e g^ ~ and go ~ A, and the total probability in the kinematical region constrained as 
go < A scales approximately as A^. However at higher g^ both the dominance of the diagonal 15 — IS 
transition in the sum in Eq. (11061) becomes weaker and the behavior of the amplitude in Eq. flBTl) derived 
for soft pions becomes questionable. It is still quite likely that with these limitations the discussed here 
approach can be used up to somewhat higher values of A: A ^ 0.8 — 0.9 GeV, than those observed in 
the pionic transitions from ip' and T' and that the Fo(980) resonance places a natural upper bound on 
the region of applicability of Eq. (!67|) . As to the contribution of higher quarkonium states in the sum 
in Eq. (ll06p . for each of these states the magnitude of this contribution relative to that of the diagonal 
transition is given by 



go 
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The transition polarizability a^^"^ should considerably decrease with n. This is supported by the 
very small experimental rate of the decay T(35') — > ttttT (a discussion of this decay in terms of the 
transitional chromo-polarizability can be found in Ref. |lU8j ). Thus, most likely, the only real effect 
of higher states up to A ~ 0.9 GeV reduces to that of the ip' resonance. This effect however can be 
accounted for in the data analysis, since for this resonances all the parameters (except for the overall 
relative sign of its contribution) in Eg. (11 08 p are known. Furthermore an observation and analysis of 
the discussed decays at higher values of qo, where the expression (!67ll is no longer valid, would be of a 
great interest for studies of the pion-pion scattering beyond the soft-pion region. 

Numerically one can estimate from Eq. (ll07p the total rate in the kinematical region constrained by 
A = 0.9 GeV as 



T{J/ij ^ 7r+7r-f+r)|go<o.9Gey ~ 10" 



TeeiJ/i^) . (109) 



2GeV- 

Given that the diagonal polarizability is likely to be larger than the transition one, it can be expected 
that for the J/tp resonance the branching ratio of the discussed decay in a useable kinematical range 
should be at the level of 10~^ which looks to be well within the reach with the expected CLEO-c data 
sample. 



5.4.3 Slow J/ip in Nuclear Medium 

Understanding the charmonium interaction with nuclear matter is important for description of the 
photo- and hadro- production of charmonium and charmed hadrons on the nuclear targets as well as 
for diagnostics of the hadronic final states in heavy-ion collisions and search for Quark Gluon Plasma. 
Such interaction has been a subject of numerous studies with a broad range of theoretical predictions. 

First perturbative QCD calculations |125lll23j predicted very small J/ip dissociation cross section by 
hadrons, on the order of few yubarn. With all the great interest to the problem of charmonium interaction 
with nucleons and nuclear matter and its practical importance, the discussion of this interaction is still 
wide open. In particular, the estimates of the strength of the interaction of J/ ip and ip' with the nucleon 
range, in terms of the scattering cross section at low energy, from a fraction of millibarn |126[ 1127] up 
to 10 mb or more |128l 11291 1130] . Recent reviews of the subject and further references can be found in 
the Refs. [ml fT32l 1133] . 

In many of these applications the most interesting energy region is usually well above the threshold, 
where the complexity of the problem becomes more confounding due to the multitude of possible 
inelastic processes contributing to charmonium scattering on nuclear matter. However the strength 
of the interaction at energy close to the threshold is also measurable [128] and its reliable estimate 
can serve as a useful reference point for analyses of the behavior of the interaction at higher energies. 
Furthermore, the J/ip and ijj' interactions at low energies are of explicit importance for high energy 
heavy ion collisions since the relative motion between the co-moving charmonium and nuclear matter is 
rather slow. Moreover the forward elastic scattering amplitude can be related to the J/ip and ijj' mass 
shift in matter predicted by a number of models [1341 11351 11361 1137] . 

The interaction of a slow charmonium with nucleons can be considered within the multipole ex- 
pansion in QCD in terms of the chromo-polarizability [127l 1138] , which parametrizes the interaction of 
charmonium with soft gluon fields inside the nucleons. For the J/ip resonance the amplitude of elastic 
scattering on a nucleon is then given by 

A{J/^ N ^ JipN) = ^ {N\ {E" ■ E'')\N) . (110) 

It can be mentioned that at low energies, below the threshold for the process J/ip + N^Ac + D, the 
only kinematically allowed inelastic reaction is J/ip + N^rjc + N. However the latter reaction involves 
the heavy quark spin-flip and should be suppressed in comparison with the elastic scattering. In order 
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to evaluate the matrix element of the gluonic operator over the nucleon in Eg. (11 101) one can make use 
of the conformal anomaly relation fIBU]) . Namely one can write 



{N\{E'^ . E^) - {B'^ . B^IN) = --{N\{G^^^f\N) = -^{Nr^\N) = ^2m% . (Ill) 

Given that the average of the operator (i?" ■ i?°) over nucleon is non-negative, one arrives at the 
inequality [138] 

{N\{E'^-E'^)\N)>^2ml . (112) 

It is expected however |138j that the chromo magnetic contribution to the nucleon mass is substantially 
smaller than the chromo-electric one, so that the actual value of the average in Eg. (11 121) should be close 
to the lower bound. 

The elastic J/ tpN scattering amplitude in the low energy limit is thus estimated as 

1571-2 

AiJ/tpN J^N) > ——a^m% (113) 

9 

and the corresponding scattering length is then given by 

ajm = ^AiJ/ijN^ J^N) "^'^ > 0.37fm (-^) • (114) 
An mj/^ + niN V2GeV / 

Accordingly, the elastic cross section at the threshold is found as 

aiJ/^pN-.J^N)=Ana'j/^^ > 17mb (^^^)' • (115) 

The positive sign of the scattering length implies attraction of the J/ip to nucleons, and the strength 
of this attraction can be expressed in terms of the binding potential for J/ip in nuclear matter with the 
number density of nucleons pN = 0.16 fm~^: 

V,, = -A(.WN^MN)£-<-2lMev{^) . (116) 

Such rather strong attraction indicates a possibility of existence of bound states of J/tp in light nuclei. 
Indeed, the condition for existence of a bound state in the approximation, where a nucleus is considered 
as being of a uniform density pat up to the sharp boundary at the radius Ra reads as 

With the minimal estimate of the binding potential in Eq. (lll6p this condition is satisfied already at 
i?yi>0.9fm, which points to a relevance of the problem of bound states to light nuclei. Although the 
criterion in Eq. (lll7p is not directly applicable for light nuclei, the resulting estimate gives credibility 
to the claims |128l 1139] that bound states of the J/ ip resonance in nuclei do exist starting from light 
nuclei. With regards to existence of a near-threshold bound or resonant state of the J/ip and a single 
nucleon, the present understanding is generally insufficient for arriving at a definite conclusion. 
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5.4.4 Interaction of slow ip' with Nucleons 

The consideration of the elastic ip'N scattering amphtude parallels that for the J/ip N process with the 
obvious replacement of by the chromo-polarizability a^i of the The latter parameter is expected 
to be larger than the 'reference' value of 2 GeV~'^, so that the numerical value of the 'minimal' binding 
energy in nuclear matter -21 MeV is very likely to be overly conservative. The shift of the mass of the 
ij)' resonance in nuclear matter can be important for consideration of the decay ip — > DD which may 
become possible due to a matter- induced shift in the mass of the D mesons P^CT| I14H [142[ 11431 11441 1145] . 

The main difference between the nuclear interactions of slow J/ and ip' is that for the latter there 
exist subthreshold scattering processes: the charm-exchange process ip'+N^Kc+D, the charmonium 
transition scattering ip'+N^J/ip+N, and generally additional channels where in the latter process 
instead of a single nucleon excited states are being produced such as A^vr, Avr, etc. The processes other 
than ip'+N^J/ip+N are beyond the scope of the present discussion. It can only be noted here that 
due to the discussed relation of the relevant gluonic matrix element to the energy-momentum tensor in 
QCD, the processes with non-diagonal transitions, such as N^Nn, should be suppressed with respect 
to the diagonal one N^N. One can also notice that similar transitions from tp' to lower charmonium 
states other than J/tp should also be suppressed in comparison with ip'^J/ip, since those other states 
cannot be produced in the second order in the leading El term of the multipole expansion. 

The process ip' + N ^ J/ip + N ioT a. slow ip' involves a momentum transfer to the nucleon ^ 
—0.82 GeV^, so that the previous, essentially static consideration is generally modified by an effect of an 
unknown form factor F{q^). For this reason the presented here estimates are somewhat approximate. 
The cross section for the discussed transitional process is found, using the known value of 
complete analogy with Eq. flllSp : 



where Pi is the cm. momentum of the initial particles. The inverse- velocity, 1/pi, behavior of the 
cross section is due the subthreshold kinematics of the process. Assuming, conservatively, that the form 
factor |F(g^)| suppresses the amplitude by not more than a factor of two, one comes to the conclusion 
that the cross section of the considered process can reach tens of millibarn at rather moderately low 
values of the initial momentum pi. 

The discussed process gives rise to a decay rate of the ip' in a nuclear medium, which can be 
evaluated |138j as 



and is likely reaching tens of MeV at the nominal average nuclear density. 

6 Charmonium above the DD threshold 

The charmonium resonances that are heavier than the D^D^ threshold at 3.73 GeV are kinematically 
allowed to decay into D meson pairs and are generally expected to be significantly broader that the 
states below the threshold. The exception from such behavior being for the resonances that might have 
mass still below the D^D*^ threshold at 3.87 GeV and the quantum numbers that forbid their decay into 
a pair of pseudoscalar mesons, such as those with unnatural spin-parity, P = (—1)"^+^, or with negative 
CP parity. Some potential models point at existence in this narrow mass range of such 'exceptional' 
resonances and both having unnatural parity. However no such states have been observed thus 
far. Instead, an already plentiful and growing suite of very interesting states is being observed at the 
£)0£)*o threshold and above, some of which almost definitely cannot be explained as simple cc states. 
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(119) 
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but rather should additionally contain light quarks and/or gluons as dynamical constituents. Thus in 
this mass range the dynamics of the heavy cc pair closely intermixes with the dynamics of charmed 
meson pairs and with general nonperturbative dynamics in QCD. 

6.1 ij{3770) 

6.1.1 General Properties 

The resonance ip{3770), or ip", with the quantum numbers J^*" = 1 can and does decay into DD 
meson pairs, which explains its relatively large total decay width[27] r['?/'(3770)] = 25.2 ±1.8 MeV. The 
e~^e~ decay width of the '?/'(3770), ree['0(377O)] = 0.247tao25 keV is about ten times smaller than that 
of the nearby Thus this resonance is considered to be dominantly a ^Di state of charmonium with 
a small admixture of ^5*1. The latter admixture enhances the e~^e~ decay rate, which otherwise would 
be very small for a pure ^Di state. The amount of mixing can be estimated from the value of Fee- 
Using a simple model with a two-state ip' ~ ^(3770) mixing Rosner[62l [63] estimates the mixing angle 
as (12 ± 2)°, which certainly agrees with the notion of the ^Di — "^Si mixing being an 0{v'^ /(?) effect 
and the estimate f^/c^ k, 0.2. It should be understood however that the particular estimated value of 
the mixing provides only an approximate guidance, not only due to its theoretical model dependence, 
but also because the experimental data, especially for '?/;(3770), are still somewhat volatile. 

In particular, the data are still not conclusive on the decay properties of V'(3770), most notably 
on the fraction of the decay rate that is not associated with the decay ip{'i77Q) DD. Namely, the 
reported by CLEQ [T46] result for the total resonance production cross section at the maximum of the 
V'(3770) peak (at E^.m. = 3773 MeV), a[e+e- ^ ^(3770)] = (6.38 ± 0.08t[i:3o) nb leaves very little if 
any room for non-DD decays, if combined with their latest measurement |147] of the DD production 
at the same maximum: o"(e^e^ DD) = (6.57 ± 0.04 ± 0.10) nb. On the other hand, the BES 
measurement of the total cross section[148j gives cr[e+e~ ■ip{3770)] = (7.25 ± 0.27 ± 0.34) nb, and 
their directly reported result |149j for the branching fractions: B['ip{S770) ~* D^D^] = (46.7±4.7±2.3)%, 
i3[V^(3770) ^ D+D-] = (36.9 ± 3.7 ± 4.2)% and B[^p{3770) non - DD] = (16.4 ± 7.3 ± 4.2)% leaves 
an ample room for non-DD decays of ?/'(3770) and their branching fraction for ?/'(3770) DD is also 
in agreement with the CLEO result for a{e^e^ DD). 

A sizable fraction of non-DD decays of ^/'(3770) would present a serious difficulty for considering it 
as a pure cc state. Indeed, its total annihilation rate should be small in comparison with such rate for 
the similarly to its leptonic width Fge. The radiative and hadronic transitions to lower charmonium 
levels are expected to be small and in fact are measured to be small: the hadronic transitions to Jijj 
all together contribute less than approximately 0.5% of the total decay rate |150l I151j . while the total 
fraction due radiative decays to 7 + r |152] likely amounts to at most about 1 %. Thus if the non-DD 
decay rate of the '?/'(3770) is measured to exceed the small rate that can be accounted for, it would 
imply an enhanced decay of this resonance into light hadrons. Such enhancement, if found, can be 
attributed [1531 1154] to a presence in the wave function of ^/'(3770) of a certain four-quark component: 
ccuu and ccdd, where the annihilation of the heavy quark pair is enhanced [155] . A mixture of the 
-0(3770) with four-quark states can be viewed as a 're- annihilation' [153] of DD meson pairs, which are 
strongly coupled to the resonance. Furthermore, within such mechanism one can expect that in the 
four-quark component an enhanced violation of the isotopic spin due to the fact that the mass difference 
between the D^D~ and D^D^ thresholds, A 9.6 MeV, is not much smaller than the excitation energy 
of the '?/'(3770) resonance above these thresholds. 

A presence of light quark-antiquark pairs in the wave function of ^/'(3770) should generally enhance [154] 
both the TTTT and rj transitions to J/ip with a larger enhancement for the latter transition, which is 
otherwise suppressed by the flavor SU(3) symmetry. Experiment ally [150j the branching fractions are 
i3[V'(3770) ^ TT+TT-J/^] = (0.189±0.020±0.020)% and i3[V^(3770) ^ r]J/^] = (0.087±0.033±0.022)%, 
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so that the ratio of the r) emission rate to that of the tt+tt^ pairs is approximately 0.5. Such ratio in- 
dicates a significant relative enhancement of the rj transition, if compared with the similar ratio, 0.1, 
for the ip' decays and given the fact that an increased phase space favors the two-pion transitions more 
than the t] emission. 

Furthermore, an isospin violation in the four-quark admixture in ?/;(3770) implies a presence of an 
isovector, / = 1, component in its wave function. Such component should then enhance the single- 
pion transition ?/'(3770) — > tt^J/i/j, and, through the ip' — iIj{3770) mixing, also provide an additional 
contribution to the amplitude of the decay tp' — n^J/ip, much needed given the previously discussed 
mismatch between the data and the theory. Starting with the needed contribution to the latter decay 
one can estimate |154j the expected rate of the former transition as B[iIj{3770) Tc^J/ip] ~ 2 x 10~^, 
which can be compared with the current experimental limit [150"] B[ip{3770) TT^J/ip] < 2.8 x 10"^ at 
90% C.L. 



6.1.2 Isospin Breaking in Production of DD pairs at iIj{3770) 

The production of the DD meson pairs at and near the ?/'(3770) resonance in e^e^ annihilation is 
essentially completely dominated by the electromagnetic current of the charmed quarks, which is a pure 
isoscalar. However the yield of the two isotopic components in the final state, D^D^ and D^D^, is 
not the same. Rather the ratio is measured |147l 1149] to be a{e~^e~ D^ D^) / a{e^e^ D^D^) = 
0.79 ± 0.01 ± 0.01 at the maximum of the '?/'(3770) peak. This is certainly not unexpected since, as 
previously mentioned, the mass difference between the charged and neutral D mesons is substantial at 
the energy of the resonance, and also the Coulomb interaction between the produced slow charged D 
mesons modifies their production cross section. Usually the effect of the mass difference in the cross 
section is estimated by the P wave kinematical factor p^ with p being the momentum of each meson in 
the cm. frame. The Coulomb effect in the limit of slow point-like particles produced by a point source 
reduces to the well known factor [1 + 7ra/(2f)], where v is the cm. velocity of either of the charged 
mesons [v ~ 0.13 for the D^D~ pairs produced at the iIj{3770) peak). A straightforward estimate of the 
product of the kinematical and the Coulomb factors gives {p+/po)^ [1 + na / {2v)] ~ 0.75 in a reasonable 
agreement with the experimental number for the charged-to-neutral yield ratio. However, it would be 
premature to conclude that the issue of this ratio is solved. Indeed, a similar estimate does not work 
at all for the B meson pair production at a similar near-threshold resonance T(4S'), where the isotopic 
mass difference for the B mesons is practically nonexistent, and the Coulomb factor amounts to about 
1.19, while the most precise data give the yield ratio very close to one [156] : 1.006 ± 0.036 ± 0.031. On 
the theoretical side, it is well understood |157l I158[ 1159] that the form factors in the production vertex 
and in the Coulomb interaction between the charged mesons generally modify the charged-to-neutral 
yield ratio. Another related effect ]T60| 1161] , that modifies both the kinematical and the Coulomb 
correction factors is the strong interaction between the mesons, which certainly is relevant since there 
is a resonance in the production channel. 

Due to small velocity of the mesons near the resonance energy one can apply the methods of the 
standard nonrelativistic quantum mechanics. With these methods the strong interaction is assumed to 
be confined to a certain radius r < a, and that in the isotopically symmetric case the (P wave) states 
of meson pairs with definite isospin, 1 = and / = 1, are characterized at r > a by the scattering 
phases 6o and 6i. Both phases behave as p^ near the threshold, while the / = phase So also makes 
a rapid variation across the isoscalar ?/^(3770) resonance. The isospin violating effects due to the mass 
difference and the Coulomb interaction (including the form factor) can be treated as being due to a 
difference 6V{r) at distances r > a in the potential for the D'^D" and D^D^ channels. Then if a 
source producing the meson pairs is localized entirely within the region of strong interaction, and the 
production amplitude has the isotopic composition A = Ao\I = 0) + Ai\I = 1), the charge-to-neutral 
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yield ratio R^/^ is given to the first order in 5V by the formula [161] 



R 



c/n 



Aq-Ai 



1 + - Im 

V 



Aq e^'^^ - Ai e^'^'' 
Ao-Ai 



1 + — 6V(r) dr 



(120) 



In the case of an essentially pure isoscalar source, relevant to the charmed meson pair production in 
e~^e~ annihilation, this general expression reduces to the following 
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1 + — 5V{r)dr 
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(121) 



so that the strong interaction effect in the discussed corrections is determined by the scattering phase 
in the isovector state 6i and does not depend on the resonant isoscalar phase Sq. 

The dependence of the effect on the parameter a is an inevitable consequence of a P wave dynamics. 
Only in the limit of vanishing phase 6i this parameter can be set equal to zero. It can be also mentioned 
that if one considers the Coulomb interaction of the charged D mesons as that of point particles, the 
effect of this interaction as well as that of the mass difference corresponds to the potential difference 
6V = A — a/r, the result of the simplified approach is recovered after integration in Eq. fll2ip down to 
a = 0: i?^/" = 1 — 3 A/ (2 vp) + 7ra/ {2v). For a nonvanishing 5i the correction depends in an essential way 
on both a and 6i )161j . Due to the dependence of this scattering phase one can expect a measurable 
variation of the ratio R^^^ with energy near the threshold. An experimental study of this variation can 
thus provide an information on the strong interaction between heavy mesons, which information would 
not be available by other means. 



6.2 X(3872) 

6.2.1 General Properties 

In the summer of 2003 the Belle Collaboration announced [162] an observation of a narrow resonance 
X(3872) produced in the decays B K X and decaying as X(3872) -n^n^ J/ip. The statistical 
significance of the new peak in the invariant mass of -k^-k^J/iI) was in excess of lOa. Shortly after 
the initial discovery the new resonance was confirmed by observation with similar significance in the 
inclusive production in pp collisions at the Tevatron [163l [TM] and by another independent observation 
in the B decays [165j . The first observed peculiar features of this resonance were its small width, 
Tx < 2.3 MeV [162j and the exceptional proximity of its mass to the threshold of D^D*^: with the 
recent improvement in the precision of the mass [166] . which placed the D^D*^ threshold at 3871. 81± 
0.36 MeV, the mass of the X(3872) corresponds to Mx - M{D^D*^) = -0.6 ± 0.6 MeV. 

The small width of X(3872) implies that its decay into DD is forbidden either by its unnatural spin- 
parity, or by negative CP parity, which, as discussed, would be possible for certain states of charmonium. 
However such states would undergo radiative transitions into 7+Xcj, which transitions were not observed 
in the experiment |162] . An analysis [168] of the angular correlations 1 1 6 7] in the process with the decay 
X 71^71^ J/ ip prefers the assignment J^'" = A similar analysis by CDF |16(I] allows either l"*""*" 
or 2^^. The latter assignment however would greatly suppress the decay of X to D^D^ir^, which is very 
close to its threshold. It is generally believed that the observation of this decay mode by Belle |170] 
and BaBar pTT] rules out the possibility of J^'-^ = . 

The positive C parity of X(3872) is in fact directly mandated by the observations [T721 1173] of the 
decay X —>■ J/tp. When combined with the existence of the discovery mode X —>■ -K^Tx'J/ip this 
implies that the pions in the latter decay have to be in a C-odd state, and thus the total isospin of the 
pion pair has to be equal to one! In particular no emission of a 7r°7r° pion pair can take place. Thus the 
X resonance definitely cannot be a pure cc system, but has to contain light quarks in its wave function. 
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Furthermore, the Belle data[172j also indicate that the decay X(3872) —>■ tt^tt^it^ J/ijj has a rate 
approximately equal to that of the decay into n~^7i~ J/ip. By the G parity a system of three pions in a 
state with a fixed C parity cannot have the same isospin as a system of two pions. Specifically, at the 
negative C parity the only possible values of the isospin are I(2tt) = 1, and /(Stt) = 0, or 2. Therefore, 
not only the isospin of X is nontrivial, but it is not definite altogether. This is also supported by the 
negative results of the search [1 74 j for charged states, which would be the isospin partners of X if it had 
1 = 1 and the isospin was a good quantum number. 



6.2.2 X(3872) as a (Dominantly) Molecular State 

The unusual properties of the X(3872) state gave rise to the suggestion [T7S1 1176^ I177i I178j that the 
wave function of this state has a significant 'molecular' component made out of mesons rather than 
out of quarks. In particular the quantum numbers J^'" = imply this component in X contains 
the state D^D*^ + D^D*^ in the S wave, which is quite natural, given the extreme proximity of the 
mass of X to the corresponding threshold. Furthermore, the threshold for the pairs of charged mesons, 
£)+£)*- _|_ £)-£)*+ is heavier by A 8 MeV, and this mass gap is large in the scale of the possible 
binding energy w for the neutral mesons. For this reason the isospin is badly broken in X(3872) and 
the wave functions of the D^D*^ + D^D*^ and D^D*^ + D^D*^ components are significantly different, 
which then explains the unusual isotopic decay properties of the X resonance. 

An existence of molecular states of loosely bound heavy hadrons was argued on general grounds long 
ago |155] and a molecular interpretation was considered[179j for explaining the properties of the then 



already known ip{AQAQ) resonance. The argument for the existence of the bound states is essentially 
quite straightforward: the strong force between hadrons containing heavy and light quarks arising due 
to the interaction between the light components does not depend on the mass of the heavy quark. 
Therefore at a sufficiently large heavy quark mass there inevitably are bound states in the channels 
where the strong interaction gives an attraction. However it had to be tested experimentally, whether 
the charmed quark is 'heavy enough' to form such states in some channels. 

It should be understood however that there is no reason to expect that only the molecular component 
is present in the wave function and that it determines all of the properties of the X(3872) boson. Rather 
one should consider the wave function in terms of a general Fock decomposition: 

tpx = aoi)Q + ^aii)i , (122) 

i 

where ijjQ is the state of the neutral D mesons 

po^*o ^ D^D*^)/^/2, while ipi refer to 'other' hadronic 
states. Due to the extreme proximity of the mass of X to the D^D*^ threshold, the ip^ part should 
be dominant at long distances. Indeed, assuming that the mass of X is below the threshold by the 
binding energy w: + m^^.o — Mx = w, the spatial extent of the tpo is determined as {rriD w)~^^'^ ~ 
5 fm (1 MeV/tf)^/^, and ^0 thus describes the 'peripheral' part of the wave function, in fact beyond the 
range of strong interaction. On the other hand, the 'other' states in the sum in the Fock decomposition 
(11221) are localized at shorter distances and constitute the 'core' of the X(3872) wave function. In 
other terms, one may think of this picture as that of a mixing in X(3872) of the molecular component 
DOj^*o _|_ /}*o/)o ^ji;]^ 'other' states, such as e.g. a 'pure' cc charmonium, which then has to be in a ^Pi 
state, also favored by the heavy quark spin selection rule [18Uj . The notion of X being a 'molecular' 
system is helpful inasmuch as the probability weight |aoP of the meson component ipQ makes a large 
portion of the total normalization. In particular the model of Ref. [181] includes S and D wave states 
of the neutral and charged charmed meson pairs as well as the channels pJ/ip and ujJ/tp, and estimates 
the weight factor of the S wave {D^D*^ + D^D*^) component as 70-80% at w ^ I MeV. 

Since the internal composition of the X(3872) can be quite different at different distances, one or 
the other part of the Fock decomposition (11221) may be important in specific processes. It appears that 
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the pionic transitions from X(3872) to J/ijj are determined by a long distance dynamics, where the 
D^D*^ + D*^D^ component dominates, so that the isospin states are mixed, and the vr^'^Tr" and t^^tt'tt^ 
transitions have approximately the same strength. The production of X however is determined by 
short distances, and proceeds through the core component, which is approximately an isospin singlet, 
as evidenced [182] by a comparable relative rate of the decays X and X K'^, while 

an exclusive contribution of the molecular D^D*^ + D*^D^ state would likely correspond to a strong 
suppression [I83l [184] of the B° ^ X decay in comparison with B^ — * X K~^. 

The mesons in the D^D*^+D*^D^ component move freely beyond the range of the strong interaction, 
where their wave function in the coordinate space is given by 

exp(-K„r) 

(pn{r)=c , (123) 

r 

where k„ is determined by the binding energy w and the reduced mass rrir ~ 966 MeV in the D^D*^ 
system as k„ = ^2 rrirW. The normalization coefficient c determines the statistical weight of the 
DO]j*o _|_ d*0£)0 component in X(3872), and its definition is correlated with that of the coefficient Qq 
in the Fock decomposition (I122p . We resolve this ambiguity in the definition by requiring that the 
coordinate wave function of the neutral meson pair be normalized to one, so that the statistical weight 
of the state {D'^D*'^ + D*^D^)/^/2 is given as |aop. If the wave function of the form (11231) is used down 
to r = 0, this requirement corresponds to c = \J k„/ (27r). 

The dominance of the D^D*^ + D*^D^ at long distances translates into a substantial isospin violation 
in the processes determined by the 'peripheral' dynamics, examples of which are apparently the observed 
decays X ti^ti^J/iIj and X -k^-k^-k^J/iIj. It is quite likely however that this isospin-breaking 
behavior is only a result of the 'accidentally' large mass difference A ^ 8 MeV between D^D*^ and 
D*0£)*o Therefore it is natural to expect that at shorter distances within the range of the strong 
interaction the isospin symmetry is restored and at those distances the wave function of X(3872) is 
dominated by / = 0. In this picture the wave function of a D^D*'' + D~~D*^ state within the region 
beyond the range of the strong interaction can be found from Eq. (11231) by requiring that at short 
distances the pairs of charged and neutral mesons combine into an / = stat^llj so that the wave 
function of the charged meson pair has the form 

, , , exp(— Kct) 

(t),{r) = c^ , (124) 



where Kc = \j2mr (A + w) k, 125 MeV. It should be noticed that both the neutral (Eq. (ll23[ )) and the 
charged (Eq. (ll24p ) meson wave function have the same normalization factor c (determined by k„) and 
differ only in the exponential power. This simple picture allows one to estimate the relative statistical 
weight of the charged and neutral D meson components in the X(3872): 



A 



(125) 



Clearly, the wave functions in Eq. (11231) and Eq. (ll24p cannot be applied at short distances in the 
region of strong interaction, where the mesons overlap with each other and cannot be considered as 
individual particles. In order to take into account this behavior an 'ultraviolet' cutoff should be intro- 
duced. One widely used method for introducing such cutoff is to consider the meson wave functions 
only down to a finite distance tq, at which distance the boundary condition of the state being that 
with / = is imposed. An alternative, somewhat more gradual cutoff, described by parameter A, can 



^The effects of the Coulomb interaction between the charged mesons are neglected here. 
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be introduced[184j by subtracting from the wave functions fll23p and (11241) an expression ce~^~^ /r. It 
should be noticed, that such regularization also results in a modification of the normalization coefficient 
c, which for the gradual cutoff takes the form 



V In A — K„ 

One can also readily see, that an introduction of any such cutoff eliminates relatively more of the charged 
meson wave function than of the neutral one, thus reducing the estimate of the relative statistical weight 
as compared to that in Eg. (11251) . so that Eq. (ll25p gives in fact the upper bound for the ratio. 

6.2.3 Peripheral Decays to D^D^tt^ and DD'j 

The 'molecular' component of the X(3872) dominating at large distances, the periphery, should give 
rise to decays to D^D^tt^ and _D°_D '^o/ [l78lll85j . The underlying processes in these decays are the decays 
of the very weakly bound D*^ meson D*^ — »• D^tt^ and D*^ —>■ D^j, and the charge-conjugate decays 
of the D*^. These decays are relatively well studied for the D* mesons, and the relevant rates can be 
deduced from the data in the Tables [27j as 

= T{D*^ D\°) = 43 ± 10 keV and = r(D*° ^ D%) = 26 ± 6 keV . (127) 

There is an interference between the amplitude of the decay of D*^ and D*^ in state with a fixed C 
parity of the initial meson pair. For the C-even X resonance the sign of the interference is positive for 
the decay into D^D^ir^ and is negative for the decay into D^D^j. The differential over the Dalitz plot 
rate of the decay X(3872) D^D^n^ can then be found as [178] 

dT{X D'D%') = |aoP ^^^f t f + '^(^"'2)1' , (128) 

where, Pq = 43 MeV is the momentum is the pion momentum in the decay of a free D*^ meson, qi and 
02 are the momenta of the and mesons in the rest frame of X, and (j){q) is the momentum-space 
wave function corresponding to that in Eg. (11231) : 

m = 4^ ■ (129) 

The significance of the interference and the total rate of the decay depend on the binding energy w: 

r(X ^ D^D%^) = \ao\^ [A{w) + B{w)] , (130) 

where A{w) describes the incoherent contribution of the decays of individual D*^ and Z)*°, and B{w) 
describes the effect of the interference between these two processes. The result of a numerical calculation 
of the terms A and B with the wave function from Eq. (ll29p is shown in FigJSl It is seen from the plot, 
that the interference between the two wave functions in Eg. (11281) significantly enhances the decay from 
a C-even state. 

The decay X(3872) DD'-f can be considered in a similar way, and the differential decay rate is 
found as 

' fl31) 



2+P -^\2-P 



(27r) 



where uo = 137 MeV stands for the photon energy in the corresponding decay D*^ D7, u is the 
energy of the photon, and p = {po — Pd)/'^ is the momentum of the D meson in the cm. frame of the 
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Figure 5: The non-coherent contribution A{w) (sohd hne) and the interference term B{w) 
(dashed) as defined in Eg. (11301) . calculated by a numerical integration in Eq. fll28p . 



final DD meson pair. An integration with the wave function with the cutoff A at short distances gives 
the expression for the total decay rate 



r(X ^ D°D°7) = r^al 



2Kn A (A + K„) / UJQ UJQ UJQ 

1 arctan h arctan — ; — 2 arctan 



UJq (A — /t„)^ V 2K,n 2 A A + 



■ (132) 



At ujq = 137 MeV and /€„ ~ 44 MeV (corresponding to the binding energy w = 1 MeV) the numerical 
value of the expression in the square braces, the interference factor, varies between 0.36 at A ^ oo and 
0.61 at A = 200 MeV. The spectrum of the photon energies is shown in FiglHl As can be expected on 
general grounds this spectrum peaks near the energy corresponding to the decay D*^ — > D^'j with the 
spread induced by the slow 'Fermi motion' of the initial meson in the loosely bound state. 

In the peripheral contribution the decay to D^D^'y is greatly suppressed in comparison with D^D^'-). 
The suppression results from three contributing factors p^5] : the small rate of the decay D*^ D^l-, 
the relatively small statistical weight of the pair of charged mesons in X(3872) and the stronger negative 
interference for the charged mesons located at shorter distances within the X resonance wave function. 
One can also argue |185j that the appearance of the pairs of charged mesons as a result of rescattering 
D^D^ D^D^ should be only a minor effect. 

A quite different final state composition should be expected for the DD'-y final state from the 
'core' component of the X(3872). In particular this component can give rise to decays through the 
intermediate '?/'(3770) resonance: X(3872) — '?/'(3770)7 — > DD'^ in which case the photon spectrum 
should have a peak at ~ 100 MeV and the isotopic composition of the final DID state should be the 
same as in the decay of '?/'(3770). 



6.2.4 e+e" ^7^(3872) as an Alternative Source of X (3872) 

The resonance X(3872) is observed experimentally only in the decays of B mesons B ^ X K and in 
inclusive production in proton - antiproton collisions at the Tevatron. Both these types of processes 
present significant challenges for precision measurements of the parameters of the resonance. In par- 
ticular, neither the total width of X(3872) is yet resolved (the current limit is Tx < 2. 3 MeV), nor its 
mass is known with a precision sufficient to determine the mass gap w from the D^D*^ threshold. A 
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Figure 6: The photon spectrum in the decay X(3872) — D'^D^'y for k„ = 44MeV (sohd) 
and for Kn = 24MeV (dashed), both at A ^ oo. The vertical scale is in arbitrary units. 




Figure 7: The production process e+e" D*^D*° 'yX{3872). 



viable alternative source of the X(3872) can be provided by the process e~^e~ — >• 7X(3872)[T86] at the 
cm. energy within few MeV of the D*^D*^ threshold, where the kinematical simplicity of the process 
would hopefully allow more detailed studies of X(3872). 

The cross section of the discussed process can be estimated by using the unitarity relation and 
considering the process, shown in FigJ7|e^e^ f)*0£)*o _^ 7X(3872), and also using the amplitude of 
the known decay D*^ —>■ D^. The amplitude for production of the D*D* pair in e~^e~ annihilation at 
a small energy above the threshold E = 2M{D*^) + W can be generically written in the form 



A(e+e- -> D*'D*') = Ao (j ■ p) {a ■ b)* + j-^A^ Pk 



a-i bk + akh - - 5ik (a ■ b) 



(133) 



where j = (e-fe) stands for the current of the incoming electron and positron, p is the momentum of 
one of the mesons {D*^ for definiteness) in the cm. frame, and Aq and A2 are the factors corresponding 
to production of the vector meson pair in the states with respectively the total spin 5 = and S = 2. 
It can be also noted that the amplitude in Eq. fll33p describes the production of mesons in the P wave. 
Another kinematically possible amplitude, the F-wave, should be small near the threshold, i.e. at 
a small W. Both Aq and A2 are generally functions of the excitation energy W. Furthermore, their 
dependence on the energy near the threshold is known to be nontrivial due to the '?/'(4040) resonance |27j. 
with possible further complications in the immediate vicinity of the threshold [1871 1188]. Neither the 
relative magnitude nor the relative phase of the amplitudes Aq and A2 is presently known, but both 
of these can be measured from angular correlations |189] . These amplitudes determine the total cross 
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section for production of D D in e^e annihilation: 

J \ m j [2tiY 2tt ^ 

(134) 

where m = M{D*^), p = \p\, and C is an overall constant related to the average value of the current 
The specific value of the latter constant is not essential for the discussed estimate, since it cancels 
in the ratio of the cross sections. The latter ratio thus can be expressed in terms of the amplitudes Aq 
and A2 and of the parameters of the decay D*^ 7 D'^, the rate and the photon energy ujq, 

aAbs(e+e-^D*°^*°^7^) _, ,2 ^jm^'^n ^2 1^0 + (9/20) 

a{e+e- ^ D*<^D*o) 2ujip \Ao\^ + \A2\^ ' ^ ' 

where u is the energy of the photon, and the form factor F is defined through the wave function (f){q) 
(Eq.([T29D as 

(p ■ k) (j) I p— — I dcos9 , (136) 




with 6 being the angle between the D* cm. momentum p and the photon momentum/c. The explicit 
form of the form factor is 



c 

pu 




ip + ^nf + < / 2 , , (p+^/2)^+A^ 

-07/2)2 + ^2 1^^ ^ 4 ^ ) {p-uj/2Y + k\ 



(137) 



with the normalization coefficient c given by Eq. fll26p . 

The 'absorptive' cross section aAbs in Eq. fll35p is (most likely) not the actual value of the cross 
section, since the amplitude of the process e^e~ — 7X can receive contribution from other mechanisms. 
Nevertheless it is instructive to examine the numerical value and the behavior with energy of this 
quantity as given by Eq. fll35p . The dependence on the cm. energy of the factor (nn/p) -^^ is shown 
in Fig|8]for two representative values of the 'molecular' binding energy in X(3872), w = 1 MeV («;„ « 
44MeV) and w = 0.3 MeV (k„ 24 MeV). This factor peaks at the energy where p ^ uj/2 ^ 70 MeV. 
The appearance of such peak is easily understood qualitatively: at p ~ k/2 the meson emerging 
from the emission of the photon in D*^ — > D^j moves slowly relative to the Z)*° and forms a loosely 
bound stately The width of the peak is clearly determined by the parameter k„. 

As is seen from the plots of FiglH] the numerical value of the factor {Kn/p)F'^ near its peak is of 
order one. Another factor in Eq. (11350 . Vquioj / {2uJq) ~ VQim/ {2uJq) ~ 1.5 x 10~^, sets the overall 
scale of the discussed cross section. The factor in Eq. fll35p depending on the presently unknown ratio 
of the (generally complex) amplitudes Aq/A2, takes values between 0.34 (at Aq/A2 ~ 0.68) and 1.31 
(at ~ 1-47), and can thus be considered as being of order one. Finally, the statistical weight 

factor |aop, as discussed, is likely to be a large fraction of one. Summarizing these numerical estimates, 
the value of the ratio in Eq. (11350 at the peak can be estimated as being of order 10~^, although the 
uncertainty is presently large. 

In absolute terms, the measured |18 71 1188j cross section cr(e+e^ £)*o£)*o-j _ 4015 MeV, i.e. 
at the energy above the D*^D*^ threshold W k. 1.6 MeV is about 0.15 nb. This cross section grows from 
the threshold as p^. With this factor taken into account the peak of the quantity (TAbs(e"''e~ 7X) 
shifts to a slightly higher value of p ~ 100 MeV, corresponding to ~ 5 MeV, where it should be 
numerically of the order of 1 pb. 

The considered mechanism of the process e^e'^ 7X describes a 'soft' production of its peripheral 
£)0£)*o _|_ D*of)0 component in radiative transitions from slow D*^l)*^ pairs. Other intermediate states 



^The same situation arises at j3 « —k/2 for the D° meson emerging from — > D^j. 
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Figure 8: The factor UnF'^/p vs. the excitation energy W above the D*^D*^ threshold at 
representative values of the binding energy w in X(3872) and the ultraviolet cutoff parameter 
K: w = IMeV, A = 200 MeV (solid), w = 1 MeV, A = 300 MeV (dashed), w = 0.3 MeV, 
A = 200 MeV (dashdot), and w = 0.3 MeV, A = 300 MeV (dotted). 



with charmed meson pairs, i.e. D^D*^ + D^D*^, DD and DD* {DD*), can potentially contribute to 
the discussed process e'^e" •yX. However, one can readily see that in either of these processes the 
charmed meson emerging after the emission of the photon is very far off the mass shell in the scale 
of Kn- Thus none of such processes can proceed due to the long-distance peripheral component of the 
X(3872) resonance, but rather is determined by the short- distance dynamics of the 'core' of X. For 
this reason such contributions, as well as other possible mechanisms related to the 'core' dynamics, 
should be smooth functions of the cm. energy on the scale of few MeV around the D*^D*^ threshold, 
where the considered absorptive amplitude experiences a significant variation. Therefore even under 
the most conservative (and quite unlikely) assumption that these mechanisms cancel the contribution 
of the latter amplitude near its maximum, such cancellation cannot take place at all energies in the 
considered range. Thus the cross section of the process e^e~ — > 7X at an energy within few MeV of 
the D*^D*^ threshold has to be at least as large as the estimates for ciAbs near its maximum i.e. of the 
order of 1 pb. The latter is a conservative estimate, since one cannot exclude that the contribution of 
those 'other' mechanisms exceeds the calculated amplitude and that the actual cross section is larger 
than (TAbs- 

6.2.5 One- and Two-Pion Transitions from X(3872) to Xcj 

Although the bulk of the data on X(3872) indicate that it is very likely to be related to dynamics of 
DOjj*o cJiaj^jxigfi meson pairs, a possibility is still being considered [1201 [121] that the observed properties 
of X(3872) can be, to an extent, mimicked by a 2^Pi state of charmonium, so that any 'molecular' 
admixture would be viewed as a secondary effect due to the coupling to the DD* states. In this picture 
the main available indicator of a significant isospin violation in X(3872), the approximately equal rate 
of the decays X(3872) pJ /''p — ^ ti^-k^ J /ip and X(3872) — > ujj /ip — > ti^ti^tx^ J /ip is explained by the 
kinematical suppression of the isospin-allowed transition X(3872) — >• uiJ/ip. 

The transitions from X(3872) to the XcJ charmonium states with emission of one or two pions, which 
can be studied in addition to the observed processes X(3872) — >• Ti'^7i~J/ip and X(3872) — > TT^7i~7i^J/ip, 
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may prove to be instrumental in further exploration of the X resonance. Such transitions may be 
accessible for experimental observation and may hold the clue to understanding the isotopic structure 
of the X(3872) and of the prominence of the four-quark component in its internal dynamics. The 
characteristics of such transitions are generally completely different between the possible charmonium 
and molecular components of the resonance X(3872). The rate of the one-pion transition relative to the 
process with two pions is sensitive to the 1 = 1 four-quark component of the X(3872), while the isoscalar 
four-quark component should give rise to relative rates of two-pion transitions to the XcJ states with 
different J, which are very likely at variance from those expected for transitions between charmonium 
levels [192]. 

The pion transitions from a charmonium excited ^Pi state are described by the multipole expansion. 
The two-pion emission is not suppressed by the isotopic symmetry and dominantly proceeds to the Xci 
state. Indeed, this is the only final state in the transitions from ^Pi where the pions can be emitted in 
the S wave, for which the pion pair creation by the gluonic fields is enhanced by the conformal anomaly 
in QCD (Eq. fp7|) ). A transition to the Xc2 resonance requires a presence of a D wave, proportional to 
a small parameter k 0.2, and is also suppressed by smaller available phase space. Numerically, the 
suppression turns out to be quite strong, about a factor of 10~^ |192j . Moreover, due to the spin-parity 
properties the Xco resonance cannot be produced in either S or D wave dipion emission, so that such 
transition arises only starting with the fourth power of pion momenta in the chiral expansion and should 
also be quite small. Furthermore, if the quarkonium matrix element (the chromo-polarizability) for the 
2^ Pi l^Pi transition is evaluated from the known amplitude of the transition —>■ ttttJ/iIj, the 
transition rate can be estimated [192] as r(2^Pi — > tt+tt" l^Pi) ~ 1.5 keV, which is most likely only a 
tiny fraction of the total width of X(3872). 

The isospin-violating single pion transitions in charmonium, 2^Pi l^Pj can also be described 

within the multipole expansion, similarly to ip' J/'ip- The ratio of the transition rates Tj to final 

states with different J is then given by 

r2 : Ti : To = 3^,^(2) : 5^,^(1) : ^ 1 : 2.70 : , (138) 

where (j) stands for the pion momentum in the corresponding process. One can notice that the Xco 
final state is not accessible in this transition too. It can be noted however that in the case of pure 
charmonium the single pion transition should be even weaker than the two-pion: in the ratio of the 
rates the unknown quarkonium matrix element cancels, and one finds [192] 

^ ^ ' ^ , ~ 0.04. 139 

A somewhat different pattern of the transition rates can be expected if the resonance X(3872) is 
dominantly a molecular state or, generically, is a four-quark state. Then the pion transitions to the 
charmonium states Xcj can be treated as a 'shake off' of the light quarks. In particular, the spin 
dependent 'heavy-light' quark interaction is proportional to the inverse power of the heavy quark mass 
mg^, so that any exchange of the polarization between the light and heavy degrees of freedom is expected 
to be suppressed. Neglecting such exchange one can arrive at the following estimate of the relative rate 
of the single-pion transitions: 

To : Ti : r2 = 4p^^,^ : Sp^^^^ : 5p^^,) ^ 2.88 : 0.97 : 1 , (140) 

so that unlike for a pure charmonium the rate of the transition to the lowest Xcj state should be 
the largest. Furthermore, due to the apparent strong isospin violation in the wave function of the 
X(3872) resonance, the single pion process should not be suppressed by the isospin, and in fact is likely 
to dominate over the two-pion transitions. For the relative strength of the latter decays to Xcj with 
different J, it can be mentioned that the transition to Xco is still suppressed in the chiral expansion by 
the spin-parity properties, while the Xc2/Xci ratio is uncertain with the kinematics obviously favoring 
the Xci final state. 
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6.2.6 X(3872) as a Virtual State 



The expected decay X(3872) —>■ D^D^tt'^ was sought for and observed |170i 1171] , however the charac- 
teristics of the observed process do not quite look like what one would expect for decay of a bound 
state. Namely, the experimental study of the B meson decays B D'^D^t^^ K |170[ 1171] and B 
D^D'^'j K |171j revealed that the invariant mass recoiling against the Kaon displays a significant en- 
hancement with a maximum at approximately 3875 MeV, which is only about 3 MeV above the D^D*^ 
threshold. The observed events can all be in fact attributed to the process B [D^D*^ + DD*^) K 
since no distinction between the D*^ mesons and their decay products was done. Moreover, the yield 
of the heavy meson pairs within the above-threshold peak is about ten times larger than that of the 
-K^-K^J/il) and n^-K^TT^J/ip channels at the peak of X(3872). It has been most recently argued [193] 
that a very plausible explanation of the observed enhancement of the D^D*^ production combined with 
the smaller observed X(3872) peak in the Ti^-K^J/ip channel is that both these phenomena are due to 
a virtual state |194l 1195] in the D^D*^ channel. In this picture the observed peak in the Ti^-K^J/ip and 
TT^ TT~ TT^ J / %p mass spectra is in fact a cusp with a sharp maximum at the D^D*^ threshold. 

The isospin properties of such near-threshold virtual state can be analyzed |196j within the ap- 
proximation of small interaction radius. Such approach is similar to the 'universal scattering length' 
approximation [TS3], and differs in including the effect of the nearby threshold for charged mesons 
D^D*~ . An interesting energy-dependent behavior of the isotopic properties arises from the mere 
fact of the mass splitting A between the two isospin-related and coupled DD* channels. In partic- 
ular it can be argued that the expected pattern of the isospin breaking is consistent with the ob- 
served relative yield of -K^-K^J/ip and h^h^ti^J/iIj at the peak which experimentally [168j corresponds 
to B{X 71+71-71^ J/ i))/B{X n+n-J/ij) = 1.0 ± 0.4 ± 0.3. Moreover, the production amplitude for 
the 1 = 1 state tt^ti~J/iI) in the considered approximation necessarily has a zero between the D^D*^ 
and D^D*~ thresholds, thus reducing the apparent width of the cusp and putting it in line with the 
experimental limit [162] F < 2.3 MeV on the width of the peak in this particular channel. 

The approximation of small interaction radius is applicable at small energy E = M{DD*) — M (D^) — 
M{D*^) for a consideration of the strong dynamics of two coupled channels D^D*^+DD*^ and D^D*^ + 
D~D*~^, which for brevity can be called n and c channels. The energy range of interest for the discussion 
of the X peak is from few MeV below the n threshold and up to the c threshold, i.e. up to ~ A = 
M[D~^D*~) — M{D^D*^) ^ 8.1 MeV. In this range the scale of the cm. momentum (real and virtual) 
in either channel is set by ■\/2m,.A ^ 127 MeV, where ~ 970 MeV is the reduced mass for the meson 
pair. One can apply in this region of soft momenta the standard picture of the strong- interaction 
scattering (see e.g. in the textbook |197j ). where the strong interaction is localized at distances r < tq 
such that Tf) a/2?tvA can be considered as a small parameter. Considering for definiteness an energy 
value between the two thresholds, < E < A, one can write the corresponding wave functions (up to 
an overall normalization constant) as 

X„(r) = sin(/c„r + 5), xdr) = ^ exp(-Kcr) , (141) 

where kn = \p2/n\E and k^. = ^2mr{A — E), S is the elastic^ scattering phase in the n channel and 
the constant ^, generally energy- dependent, describes the relative normalization and phase of the wave 
function for the two channels. 

The wave functions (11411) should be matched at r ^ tq to the solution of the 'inner' problem, i.e. 
that in the region of the strong interaction. In the limit of small tq all the complexity of the 'inner' 
problem reduces to only two parameters. Namely, in the region of the strong interaction the n and c 
channels are not independent and get mixed. Due to the isotopic symmetry of the strong interaction 

^In this consideration the smaU inelasticity due to the tt+tt" J/V' and tt+tt^tt^ J/?/' channels is neglected and will be 
included later. Also the small width of the D* mesons is entirely neglected. The effects of the latter width are considered 
in the most recent papers (1951 1199] . 



53 



the independent are the channels with definite isospin, 7 = and / = 1, corresponding to the functions 
Xo = Xn + Xc and Xi = Xn — Xc-, and the matching parameters are the logarithmic derivatives — and 
— Ki of these functions at r = tq. Using the assumption of small the matching condition for the 
functions from Eq. (ll4ip can be shifted to r = 0, so that one can write the resulting matching equations 
as 

fc„COs5-^Kc fcnCOS^ + ^Kc 

■ X I c = ' — T~ = • 

sm + 4 sm — 4 

These equations determine both the scattering phase 6 and the constant ^ as 

cot 5 = -^ (143) 

with 

«efr = ^ , (144) 



and 



"° sin5. (145) 



The nonrelativistic scattering amplitude in the n channel is therefore given by [197] 

1 



5 



and the scattering length a is thus found from the E = limit of this expression as 



(146) 



1 



Ko + V2m^A ^^^^^ 



E=o 2KqKi - {kq + Ki)y/2m,A 



The whole approach is applicable if the scattering length is large in the scale of strong interaction. 
A large positive value of a implies an existence of a shallow bound state, while a large negative a 
corresponds to the situation with a virtual state |197] . According to the estimates of Ref . |193] the 
required by the data scattering length in the problem considered is — (3 -i- 4) fm, corresponding to a 
negative and quite small indeed parameter hies{E = 0) ~ (50 60) MeV. 

The physical picture, consistent with a small ^eff ) and which could be argued on general grounds [155l| . 
is that an attraction in the 7 = channel is strong enough to provide a small value of kq, while the 
interaction in the 7 = 1 channel is either a weak attraction or, more likely, a repulsion. In both cases the 
absolute value of ki is large, i.e. of a normal strong interaction scale, with the sign being respectively 
negative or positive. Another, purely phenomenological, argument in favor of large \ki\ is that no 
peculiar near-threshold behavior is observed in the production of the 7 = 1 charged states, e.g. D^D*^ . 
At large the expression fll44p simplifies and takes the approximate form 

Keff ~ 2ko - . (148) 

Using this approximation, one can readily see that in order for Kes{E = 0) to be negative and small, 
the parameter kq has to be positive and quite small: 



Ko < ^ 63 MeV. (149) 

It is interesting to note that in the discussed picture the interaction in the 7 = state is strong 
enough by itself to produce a shallow bound state in the limit of exact isospin symmetry, i.e. at A — 0. 
In reality the isospin breaking by the mass difference between the charged and neutral charmed mesons 
turns out to be sufficiently significant to deform the bound state into a virtual one, i.e. to shift the 
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pole of the scattering amplitude from the first sheet to the second sheet of the Riemann surface for the 
amplitude complex function of the energy E. 

The inelasticity in the n and c channels, related to decays to the observed final states tt+tt^J/-?/' 
{pj/ip), n^n^n^J/ip (ujj/ip), ■yJ/ip and probably other, appears to be reasonably small, as one can infer 
from the observed |170l 1171] dominance of the D^^D*^ production in the threshold region, and can be 
parametrized by a small imaginary shift 27 of the denominator of the scattering amplitude in Eq. fll46p : 

(150) 



If one further assumes |18Ul 1193] that the 'seed' decay B — * XK is a short-distance process, one would 
find that the yield in each final channel coupled to X is proportional to that channel's contribution to 
the unitary cut of the amplitude F. This implies in particular that 

^[5 ^ (D°D*° + 5/^*°) i^] : B{B^uJ/^K) : 13{B ^ pJ/^K) = kn\F\^ : Ju\F\^ ■ lp\FW 

(151) 

where the specific expression for |Fp depends on the value of the energy E relative to the n and c 
thresholds, as given by Eg. (11 50 p an its analytical continuation across the thresholds. Besides the energy 
dependence of the overall factor the heavy meson channel contains the phase space factor while 
for the ujJ/ip and pJ/ip yields an additional dependence on the energy arises from the factors 7^^ and 
7p- 

A certain variation of the width parameter 7^ for the tt+tt ix^J/ip channel in the discussed range of 
energy is of a well known kinematical origin. Indeed, the central value of the mass of the uj resonance 
puts the threshold for the channel uJ/if) at 3878.5 MeV, which corresponds to -E ~ 6.7 MeV in our 
conventions, i.e. squarely between the n and c thresholds. Any production of the -k^-k^-k^J/iI) states at 
smaller invariant mass is a sub-threshold process, possible due to the width F^^ of the uj resonance. In 
other words, the energy dependence of the width factor 7^ can be estimated as 

^^ = \A^W^^ , (152) 

where is the amplitude factor for the coupling to the ujJ /i}) channel and is the effective momentum 
of UJ at the invariant mass M calculated as 

QesiM)= \q{m)\— 153 

with the cm. momentum |(f(m)| found in the standard way: 

^[{M - mji^Y - m2] [(M + mji^Y - 

m'^)\ = - ym ■ ^ ^ 

The lower limit in the integral in Eg. (11531) can be chosen anywhere sufficiently below — T^, since 
the Breit-Wigner curve in the integrand rapidly falls off away from the resonance. 

Numerically, the effective momentum ([^^ can be estimated as varying from approximately 20 MeV 
to 50 MeV between the n and c thresholds, i.e. when E changes from = to = A. In the pJ/V' 
channel the expected energy behavior of the yield is guite different. If one writes the corresponding 
width factor 7p similarly to Eg. (11521) as 

7p=|A,pgi^ , (155) 

the effective momentum can be estimated as varying only slightly due to the large width of the p 
resonance: ~ (125 -h 135) MeV as the energy changes between = and E = IS.. 
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Figure 9: The expected shape (in arbitrary units) of the virtual state peak in the yield of 
Ti^-K^ J/ip (solid) and -k^ -k^ -k^ J / ip (dashed) channels. 



The amplitudes A^^ and Ap on the other hand should display a noticeably different behavior in the 
energy range of interest, due to the rapid (and different from each other) variation of the / = and 
1 = 1 scattering amplitudes of the meson-meson scattering. Namely assuming that the amplitudes for 
production of the / = and / = 1 at short distances are described by constant (in the considered energy 
range) factors and $p and that the effective strong production radius is i?, one can arrive [T96]|, after 
taking into account the rescattering of mesons, at the following estimate for the ratio of the amplitudes 

^='^^{l-n.R)^. (156) 

Since the situation where the X peak is a virtual state corresponds to a small positive kq satisfying the 
condition (11491) . the amplitude Ap described by Eq. (ll56p should necessarily change sign between the n 
threshold, where Kc = y/2'mrA, and the c threshold, where Kc = 0. 

The expected difference in the shape of the cusp in the pJ/ ip and uJ/ ip channels is illustrated in Fig|9l 
In these plots the parameters of the virtual state correspond to the scattering length a = —(4 + 0.5 i) fm, 
which is close to the possible fit values of the scattering length found in Ref. |193j . In the limit of large 
Ki this value of a translates into kq ^ 38 MeV and 7 6 MeV. One can see from Figl9] that due to the 
discussed zero of the amplitude, the peak in the ti^7i~J/i/j channel is expected to be quite narrow in 
agreement with the experimental limit on the width of X(3872). The plots in Figj9]are normalized to the 
same total yield in each channel over the shown energy range in order to approximate the experimentally 
observed relative yield. Such normalization corresponds to setting 



1 - KiR 



175 MeV 



which value does not appear to be abnormal, even though at present we have no means of independently 
estimating this quantity. 

Summarizing the situation with the X(3872) resonance it can be stated that as of the time of this 
writing, full four years after its discovery and even though this resonance is listed among established 
particles in the Tables [27] . its real status is still lively debated in the literature. It is clear that this 
peak is very closely related to the near-threshold dynamics of the DD* meson pairs, but it is still 
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unclear whether this is a shallow bound molecular state, or a virtual state, or something else. It may 
well be that a further pursuit of the puzzles arising in connection with the X(3872) peak may provide 
important clues to understanding multi-quark dynamics. 

6.3 Higher J^^ = 1 Resonances 

The existence of reasonably broad resonances above the open charm threshold has been fully expected on 
the basis of potential models[15j. What has not been fully expected however are the peculiar properties 
of individual resonances, different for different states, and never failing to bring unexpected surprises. 

6.3.1 The Vicinity o/^(4040) 

An unusual behavior of the cross section of e~^e~ annihilation into charmed meson pairs in the energy 
region, which is now associated with the resonance ip{4040), was first pointed out in Ref. |179j . Namely, it 
was noticed that the production of the D*D* is greatly enhanced relative to DD and DD* in comparison 
with a simple spin-model estimate [179], [15]. The cross section in each channel is proportional to the P 
wave factor with p being the cm. momentum in the corresponding final state. The model predicted 
the ratio 1:3:7 of the extra factors on top of the p^ in the cross section for production of respectively DD, 
DD* + DD* and D*D* . The then existing data however indicated that the cross section for production 
of the pairs of vector mesons D*D* near their threshold was enhanced by a factor of tens to hundreds in 
comparison with this estimate. This observation gave the reason for the suggestion [TTQ] that ^/^(4040) is 
in fact a molecular state made of the vector mesons. However a more conventional potential model [15] 
could well accommodate the ip{AQ4Q) peak as a dominantly S^S*! state, and it was also argued (200j that 
the suppression of the decay modes for this resonance with one or two pseudoscalar mesons is due to 
an 'accidental' kinematical zero of the overlap integrals in a model of such decay. 

A significantly more detailed, than previously available, data on the production of charmed meson 
pairs have been accumulated recently ^lST] 1188^ I2U1] and the new data also include the cross section for 
production of pairs of strange charmed mesons DsDs- The new data still indicate an unusually strong 
production of the vector mesons, D*D*, near their threshold and also reveal intricate features of the 
behavior of the cross section in other channels as the energy sweeps across the D*D* threshold(s) and 
the ip{4040) resonance, which definitely points at a strong coupling between the channels. This behavior 
is more illustrative in terms of dimensionless rate coefficients Ri defined as follows [202] 

D,Ds)=Ms)vl^R2 , (157) 

D*D*) = aois)7iv^o + vl)R, , 

with (To = 7rQ;^/(3s). Here vd, vn^, Vq and f+ stand for the cm. velocities of each of the mesons in 
respectively the channels DD, DgDg, D*^D*^ and D*^D*~, while for the channel DD* + D*D with 
mesons of unequal mass the velocity factor is replaced by {2p/ y/s) with p being the cm. momentum. The 
values of Ri calculated from a preliminary version of the data |187] are shown as data points in FigJTOl 
The extra factors 1, 3 and 7 in Eq. (11571) for respectively the pseudoscalar-pseudoscalar, pseudoscalar- 
vector and the vector-vector channels correspond to the ratio of the corresponding production cross 
section in the simplest model[T79l [15] of independent quark spins, so that the inequality between Ri, 
i?3 and i?4 also illustrates a conspicuous deviation from this model. In particular the very large values 
of i?4 describe the unusually strong enhancement of the vector-vector channel. 

The curves in the plots of Fig{TO] correspond to a fit |202j with one resonance and a non-resonant 
background with coupling among the channels and with a proper treatment of the onset of the inelas- 
ticity at the two thresholds for the vector meson pairs, D*^D*^ at 4013 MeV and D*"^D*~ at 4020 MeV. 



a(e+e- ^ DD) = (To(s) 2v]^Ri , 
a{e-^e- ^ DD* + D*D) =aQ{s)Q (^^^ Rs , 
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Figure 10: The plots of the rate coefficients R corresponding to the fit with excluded data 
points at energy 4010 and 4015 MeV for the DD channel and at 4015 MeV for the D*D* 
channel. The excluded points are shown by filled circles. 



A fit to the full set of data between 3970 MeV and 4060 MeV turns out to be impossible with an ac- 
ceptable (the best fit corresponds to /NDF = 17.6/8. The offending are the data points at 4010 
and 4015 MeV for the DD channel and, to some extent, the data point at 4015 MeV for the D*D* 
production. If those points are removed, a fit to the data can be found with 'x^ /NDF = 3.0/5 and the 
resulting central values of the resonance nominal mass, width and Fee are M = 4019 MeV, F = 65 MeV 
and Fee = 1.7 keV. A more detailed analysis [203] of the final data including the radiative corrections 
resulted in the values M = 4013 ± 4 MeV, F = 66 ± 8 MeV and Fee = 1-9 ± 0.7 keV. 

The peculiar behavior of the DD points at 4010 - 4015 MeV may indicate a presence of another 
narrow resonance at this energy at or between the ]J*^l)*^ and D*^D*^ thresholds, which has a quite 
small coupling to e~^e~, corresponding to Fee in the range of a few tenths of eV|202]. Such resonance 
could then be explained as a P wave D*D* molecule which is expected |155] to have the e"'"e~ width in 
the same ballpark. 

6.3.2 ^(4170) 

The main surprise of the peak ip{4:170) (also labeled as ■?/'(4160) in the Tables [27]) is that it corresponds 
to the absolute maximum of the cross section for production in e'^e~ annihilation of strange charmed 
meson pairs DsD* + DgD*. The cross section at the maximum reaches fl87l 1188] nearly 1 nb, which is a 
very large value for an exclusive state with hidden strangeness and charm. This property of the peak 
is intensively used in the experimental studies of the Dg mesons, but it has no theoretical explanation. 

6.3.3 r(4260) 

The total cross section for charm production in e^e~ annihilation has a well known [27] deep minimum 
around the cm. energy 4260 MeV. It is at about the same energy where the BaBar experiment, while 
performing a survey of the region around 4 GeV in the radiation return process e~^e~ — 7 X, found [20^ 
a peak in the cross section for the exclusive final state Tt^-K^J/ip. The existence of the peak was 
further con&med by Belle |205] and by CLEQ [206] using the same radiative return method and also by 
CLEO [207j by a direct scan of the e"'"e~ annihilation in the peak region. The peak at the same invariant 
mass is also indicated by the data [182] on the decays B Tr~^7T~ J/ipK. The values for the mass and 
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the width measured by independent experiments are in a statistical agreement with each other, and the 
average values of these parameters according to the updated Tables ^27j are M(Y) = 42641*^12 MeV and 
r(y) = SStlj MeV. 

The only so far observed decays channels of y(4260) are tt^ti^J/iIj, -K^-K^J/ip (in a fair agreement 
with the isospin of Y being equal to zero) and K'^K^ J/ip. Notably, the data show no peak in any 
channels with D meson pairs. In particular the most stringent upper limit is found [208] for the DD 
channel: r{Y — > DD)/r{Y — > ir^ 7r~ J / ip) < 1.0 at 90% confidence level. Such behavior is next 
to impossible to explain by considering the resonance as a charmonium state even though the mass 
of Y is close to the expectation for the 45" state [209] . since lower J^'" = 1 resonances ip{3770), 
V^(4040) decay practically exclusively to D meson pairs. A whole spectrum of interpretations of the 
state F(4260) has been suggested |210j including cc + glue hybrids [2m 1212] and a tetraquark cscs 
state assignment [213j . which does not look very promising, given that the ratio of the decay rates [206] 
r(y — i> K~^K~ J/ip)/T(Y — > 7i~^Tr~ J/ip) ~ 0.15 does not show any enhanced presence of strangeness 
within r(4260). 

6.3.4 = 1— Peaks m n+n-ip' 

Subsequent studies of the radiative return events have lead to an observation by BaBar[214j of a 
'broad structure' in the final state TT^n'tp'. A single resonance fit to the data yielded the mass of 
M = 4324 ± 24 MeV and the width of T = 172 ± 33 MeV. A further investigation of this final channel 
resulted in an observation by Belle [2 15] of a peak with the mass of 4361 ± 9 ± 9 MeV and the width 
of 74 ± 15 ± 10 MeV, possibly compatible with the structure observed by BaBar, and an additional 
narrower peak at 4664 ± 11 ± 5 MeV with the width of 48 ± 15 ± 3 MeV. 

6.3.5 Z(4430) and Remarks on New States 

Most recently the Belle experiment presented data [2 16] on observation of a peak Z(4430) in the charged 
system n^ip' emerging from the decays B Kn'^ip'. The parameters of the peak are M = 4433 ± 4 ± 
2 MeV, r = 45t\l -11 

MeV. The statistical significance of the observation corresponds to 6.5(j, however 
in view of an utmost importance of the observed peak an additional confirmation is eagerly awaited. 
Unlike the previously discussed electrically neutral states, which all at least had a chance of being a 
pure charmonium, this one is charged and has isospin 1 = 1 and clearly cannot be a pure cc but has to 
contain light quarks in addition to the cc pair. Naturally, a variety of interpretations of Z(4430) has 
been suggested: a threshold peak |217l 1218] or a resonance [2 19] a loosely bound molecular state [220] in 
the D*Z)i(2420) meson system, a radially excited tetraquark |221j . a QCD-string based mo del [2^ . and 
a baryonium state [223] . In either case, whether or not the peak Z(4430) is related to the D*Z)i(2420) 
meson system, its 'affinity' to the particular decay channel nip', rather than the multitude of available 
channels with D mesons presents a very intriguing riddle. 

Clearly, a similar riddle also relates to the states y(4260), and the 'structures' in the n'^TT'ip' channel 
at 4.32 and 4.66 GeV. Perhaps, the simplest explanation of such unusual prominence in the decays of 
each state of a particular charmonium level accompanied by light hadrons would be a picture, where 
a charmonium state, e.g. J/ip, or ip' is 'stuck' in a light hadronic state. In a sense, such picture can 
be viewed as that of a bound state of a relatively compact charmonium inside a light hadron having a 
larger spatial size. This possibility in fact returns us to the previous discussion of existence of bound 
states of J/ip and/or tp' in light nucleons. It may well be that some of the recently found high mass 
states are in fact such bound states in mesonic rather than baryonic matter. If this indeed is the case, 
one can naturally expect an existence of similar baryonic states, e.g. bound states of either J/ip or ip' 
'inside' a proton, or even 'inside' a deuteron. Needless to mention that an observation of such baryonic 
states would be of an immense interest. 
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7 Summary 



The potential models of interaction within charmonium generally agree with the lower-mass part of 
the observed spectrum of the resonances, where two long-missing states, the ^Pi he and the 2^5*0 r)'^, 
have eventually been located. In the mass region at and above the threshold for charmed mesons their 
dynamics apparently plays an important role in determining the spectrum of states. At lower masses 
some fine effects of the interaction between the quark and the antiquark still remain unsolved, such as 
the behavior of the spin-dependent forces, describing the splittings of the XcJ states, and the spin-spin 
interaction giving rise to the small splitting between the c.o.g. of the XcJ resonances and the spin-singlet 
he- 

Less model-dependent and based on the underlying theory of QCD spectral methods for analyzing 
the charmonium states, either analytical, or numerical lattice simulations, reproduce reasonably well the 
properties of the lowest state in each channel. However, these methods are intrinsically less appropriate 
for handling radially excited states, and they also eventually run into limitations discussed in the Section 
2.3.2. 

The overall picture of the hidden charm decay through strong or electromagnetic annihilation is in 
agreement with the data, still leaving us with a number of puzzles in some particular cases, by solving 
which we might gain some new insights into hadron dynamics. Among such puzzles are the larger than 
expected total widths of r)c and Xco and the non-similarity of some exclusive decay channels for J/ ip and 
It'' - a violation of the '12% rule'. These yet to be explained properties of the well known charmonium 
states very likely indicate a presence of substantial nonperturbative effects in the annihilation processes. 

The radiative transitions between charmonium levels arc in line with generic considerations based 
on nonrelativistic quantum mechanics and also with the results of specific potential models. A notable 
exception is the Ml transition J/ip ^ 777c, which is significantly weaker than any theoretical estimates. 
Coupled with the large total width of rjc this may signal a mixing of the rjc with light J^^ = degrees 
of freedom. 

The hadronic transitions between charmonium resonances offer a window into the interaction of 
the charmonium states with soft gluon fields, and also into the details of the conversion of soft gluons 
to light mesons. The former interaction is described by the multipole expansion in QCD, while the 
latter conversion is tractable due to the chiral algebra and low-energy theorems in QCD based on the 
conformal and axial anomalies. Such approach turned out to be successful in describing the transitions 
i/j' TrnJ/ip and ip' riJ/ip and in predicting finer details such as the small D wave in the former 
transition. Due to this success it can be believed that the recently established discrepancy by a factor 
of about 1.5 between the observed and the theoretical rate of the decay ^' 7r°J/^ is not a result 
of a failure of the multipole expansion and the low-energy theorems, but rather very likely signals a 
presence of a small four-quark isovector component in the ijj' resonance. 

The knowledge of the strength of the interaction of charmonium with soft gluon fields, the chromo- 
polarizability, can be applied to considering the behavior of charmonium inside a hadronic media, such 
as the interaction of slow charmonium with nucleons. The estimates of the chromo-polarizability from 
the rate of the decay ip' — > 7nrJ/ip indicate that slow charmonium interacts quite strongly with light 
hadrons. This interaction may well result in existence of bound states of charmonium in light nuclei or 
with ordinary hadronic resonances. 

The unusual new states at and above charmed meson thresholds keep mushrooming in the most 
recent experimental data suggesting that this mass region is a real playground for much of the exotics 
that has been previously speculated about: tetraquarks, hybrids, molecules, and mixtures of those. The 
peak at X(3872) that started it all still offers intriguing puzzles, and further studies are required in 
order to conclusively assess its internal structure and even its very status as of a resonance. The newly 
found peaks in the mrJ/ip and TiTTip' invariant mass spectra and especially the latest peak Z(4430) in 
the TT^ip' channel lead us further down the path of acceptance of multiquark hadronic states being as 
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commonplace, as are nuclei if viewed as multiquark systems. 

Many of the properties of charmonium, starting with a very small width of the J /ip resonance, came 
as a great surprise when first observed. Understanding these properties and in some cases making their 
description routine, has greatly advanced the knowledge of the dynamics of quarks and gluons. Now, 
thirty three years after the discovery of J/V', the studies of charmonium keep bringing new surprises 
and puzzles, solving which will hopefully result in further advances. 
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